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ABSTRACT 

Free standing pyrolytic graphite (FG) thrust chambers have been 
evaluated for application to high energy upper stage and attitude control 
liquid rocket engines. This Phase I Analysis Report presents the results 
of stress analyses, heat transfer analyses, and a vendor product survey and 
evaluation. The results of PG thrust chamber firings with F 2 /H 2 # Fo/fiAlOlU, 
and N 2 O 4 /O .5 N 2 H 4 -0.5 UDMH propellants at thrust levels of 100, 1000, and 
5000 pounds will be reported in a subsequent technical report. 

The analytical stress calculations which included the effects of 
PG anisotropy, creep, and growth showed the structural feasibility of fabri¬ 
cating and firing free standing PG thrust chambers up to 5000 pounds thrust. 
Operation at chamber pressures up to 500 psia appears feasible. However, the 
effects of growth and creep at high operating temperatures may be critical 
and experimental data are required to establish valid design criteria. 

The heat transfer studies investigated the transient and steady 
state temperature distributions within the PG chambers for each of the pro¬ 
pellants considered. The high inside surface temperatures encountered in 
the full scale thrust chamber using fluorine propellants create a problem 
of heat soak back in the injector after the run which must be resolved in 
future designs. 

Five pyrolytic graphite vendors supplied samples of their pyrolytic 
graphite products in the form of free standing, thin walled tubes 2 inches 
in diameter and 4 inches long. A total of 78 tubes was purchased. These 
were evaluated by burst testing in order to evaluate their structural capa¬ 
bility for free standing pyrolytic graphite rocket thrust chambers. The 
pyrolytic graphite fabrication parameters which were evaluated included 
mandrel material, mandrel surface finish, substrate nucleated and regenerative 
microstructures, controlled delaminations, boron alloying, and fiber rein¬ 
forcement. 


The critical effect of nodule size on maximum tube strength was 
confirmed. Maximum tube strengths were achieved with boron alloying and with 
regenerative microstructure. A burst pressure of 1000 peig appears to be a 
reasonable target for currently available pyrolytic graphite at room temper¬ 
ature. 


The results of continuing manufacturing studies by the vendors and 
the evaluation of test firings to be made during the experimental phase of 
this program will provide a better definition of the full potential of pyro¬ 
lytic materials as a thrust chamber material. 
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I. INTRODUCTION 


A. Background 

Future missions in space using high energy propulsion systems will 
require versatile, reliable rocket engines. These requirements pose the 
challenge of developing thrust chamber materials that can cope with the 
higher temperatures and higher heat fluxes of the high energy propellants. 

Also, space maneuvers calling for throttling and multiple starts make it 
desirable to develop lightweight thrust chamber designs which are insensitive 
to duty cycle. 

One of the most unusual structural materials to become available 
for possible application to such lightweight thrust chambers is free stand¬ 
ing pyrolytic graphite (FG). This unique form of carbon, resulting from 
the pyrolysis of a hydrocarbon gas in a high temperature vacuum furnace, is 
formed by deposition on a suitable substrate or mandrel as an impervious 
layer of pure graphite which approaches theoretical density and possesses 
unusual anisotropic physical and mechanical properties. 

By using suitable starting materials in the deposition process, 
alloys and composites with the pyrolytic graphite may also be formed with 
a wide range of properties. The main limitation on the use of pure FG 
with the storable hypergolic liquid propellants has been excessive oxidation. 
However, new pyrolytic graphite alloys currently being developed possess 
improved oxidation resistance as well as increased strength. It appears 
that pure PG is the most chemically resistant material available for the 
fluorine propellant environment. 

A free standing FG thrust chamber is formed by pyrolysis of hydro¬ 
carbon gas, with subsequent deposition of the carbon on the inside surface 
of a female mandrel. The carbon deposit is later separated from the mandrel, 
resulting in a free standing shell. The deposition, or a-b plane (as shown 
in Figure l), is parallel to the mandrel surface and is a plane in which 
high tensile strength and thermal conductivity exist in all directions. 

The shear strength in the a-b plane is very low. All directions within 
the a-b plane are usually referred to as the a-direction. In the c-direction, 
normal to the deposition plane, the tensile strength and thermal conductivity 
are very low. The elastic properties also vary greatly in different dir¬ 
ections. For example, the coefficient of thermal expansion is about fifteen 
times as large in the c-direction as it is in the a-direction. As a con¬ 
sequence of the highly anisotropic nature of FG, the thermal stresses in a 
free standing PG chamber are much more complex than those encountered with 
more conventional materials which are usually isotropic. 

The structured stresses in a free standing FG chamber are fortunately 
of greatest magnitude in the same direction as the high a-direction tensile 
strength, whereas lower tensile stresses exist in the weaker c-direction. 
However, the anistropy which is responsible for the unusual properties of 
PG also creates some special structural and heat transfer problems which 
are discussed below. 
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B. Structural Problems 

A free standing PG thrust chamber is a thin vailed axisymmetric 
shell vith regions of complex curvature near the throat and perhaps also 
near the injector-chamber Joint. 

The effects of complex curvature have not been included in the 
present work, except for some approximations of residual stresses near the 
throat. Instead, the stress analysis has been made for a cylindrical shell, 
vhlch is a good approximation for most portions of the thrust chamber. It 
is expected that the cylindrical shell approximation is also good for cir¬ 
cumferential stresses at the throat, but not for axial stresses at the 
throat, which are highly sensitive to the actual shape of the throat region. 

The structural features of free standing PG thrust chambers are 
discussed below. 

1. Residual Stress 


If the temperature of an isotropic cylinder is changed from 
one uniform temperature to another uniform temperature, the stresses in 
the cylinder do not change, if the material is perfectly elastic. However, 
an anisotropic cylinder ejperiences changes in its stress distribution when 
it is taken from one uniform temperature to another. The stresses existing 
in the PG after it is cooled from the deposition temperature to room temp¬ 
erature cure termed residual stresses. 

The change of circumferential stress on the inside and out¬ 
side surfaces of a thin walled cylindrical shell of free standing PG during 
a uniform temperature change can be calculated from the following equation: 


a 


E (or - 0 ) 
a c a 


(1 - 


12 



( 1 ) 


Where 

T = Temperature change 

E = Modulus of elasticity in the a-direction 
a 

c* c = Coefficient of thermal expansion in the c-direction 
a = Coefficient of thermal expansion in the a-direction 
t = Wall thickness 
r^ * Inside radius 

v 12 « Poisson's ratio in the deposition plane 
a = Circumferential (hoop) stress 
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A theoretical estimate of circumferential residual stresses 
can be calculated from Equation (l) by assuming that the pyrolytic graphite 
is stress-free at the deposition temperature, say kXX)*F. If room tempera¬ 
ture is assumed to be 0 # F, the residual stresses vill then equal the change 
in stress for a 4000®F temperature change as calculated by Equation (l). 

The theoretical circumferential residual stresses will have their maximum 
tensile value at the inside surface of the pyrolytic graphite and will have 
their maximum compressive value at the outside surface. Residual stresses 
also arise in the radial and axial directions. 

Probably the mo 6 t important characteristic of residual stresses 
in free standing pyrolytic graphite is in the fact, as seen from Equation (l), 
that the magnitude of the residual hoop stresses increases in direct pro¬ 
portion to t/r, the ratio of wall thickness to radius. Radial stresses are 
also proportional to t/r. As a result, the value of the thickness to radius 
ratio is a critical factor in design optimization and fabrication. It has 
been found that parts with residual hoop stresses much less than the a-direction 
strength have been impossible to produce, due to many factors such as inter¬ 
action between the mandrel and the pyrolytic graphite before separation. 

A general rule of thumb is that free standing PG parts cannot 
be fabricated without cracks or delaminations if the ratio of wall thickness 
to radius (t/r) exceeds l/lO. On the other hand, the residual stresses do 
not theoretically restrict the size of free standing parts which can be made. 
Axisymmetric shapes with diameters as large as 3^ inches can be produced in 
existing furnaces. 

2. Limited Chamber Wall Thickness 

The limitation on the maximum t/r value leads to the important 
conclusion that the maximum feasible wall thickness increases with the size 
of the thrust chamber. For example, a chamber with a 1-inch throat diameter 
is limited to a wall thickness at the throat of about 0.050 inch, whereas a 
chamber with a throat diameter of 10 inches can use a wall thickness at the 
throat of about l /2 inch. 

3 . As-deposited Residual Stresses 

Several techniques have been developed for measuring residual 
stresses in free standing ' D G. It has been found that the measured residual 
stresses are considerably different than the theoretical stresses predicted 
from Equation (l) using a stress-free temperature of hOOO* F. In fact, there 
i£ no stress-free tenperature which would lead to theoretical residual 
stresses consistent with the measured stresses. The term "as-deposited" 
refers to the residual stresses and the type of pyrolytic graphite actually 
produced by the deposition process. 

Several theoretical studies of as-deposited residual stresses 
in PG are reported in References 1 and 2. These analyses are complex, and 
their validity has not been conclusively proven. A second approach to con¬ 
sideration of as-deposited residual stresses is to rely on experimental 
data. Unfortunately, such data are available over a limited range of t/r 
values. However, the present study has used extrapolations of these data, 
as discussed in Section V. 
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4. Proof Stresses 


At motor ignition, chamber pressure will apply pressurization 
stresses suddenly to the free standing PG structure. The stresses resulting 
from superposition of the as-deposited residual stresses with the pressuriza¬ 
tion stresses, while the PG is still at room temperature, are defined as 
"proof stresses". 

5. Growth 


When pyrolytic graphite is exposed to elevated temper atures, 
in the absence of applied stress, it experiences a permanent change in dim¬ 
ensions with elongation in the a-direction and contraction in the c-direction. 
This phenomenon in PG is usually referred to as "growth", and is due to 
changes in the crystal lattice dimensions. The rate and total amount of 
growth increase with temperature, and growth is quite rapid at temperatures 
near 6000 # F (Reference 3) • 

Growth of PG in self-restrained shapes such as cylinders or 
thrust chambers has an effect on the stress distribution which is somewhat 
analogous to the effect of anisotropic thermal stresses. Strains are induced 
in the PG, with resultant changes in the stress distribution. This phenomenon 
might increase or decrease the maximum stresses in the PG, depending cn the 
temperatures, firing time, and duty cycle. 

Growth is a permanent deformation, and therefore may have 
important effects on the strength of PG either during a firing cycle, or 
even after cooldown, at which time the residual stresses will be changed 
from the original as-deposited values. An analysis of growth stresses is 
given in Reference 4. 

The effect of growth on the strength of PG is unknown. In fact, 
the strength of as-deposited PG shapes is still under continuous reevaluation. 
Therefore, it is necessary to compare the stresses after growth with the 
strength of as-deposited PG. 

6. Creep 

Creep describes the tendency of materials to permanently deform 
during exposure to stress and temperature. Although the distinction between 
creep and plastic deformation is hard to define in practice, the deformation 
attributable to creep is ideally that related to the duration of exposure 
to stress and temperature, while plastic deformation can be idealized as 
instantaneous permanent deformation. Similarly, the distinction between 
creep and growth, in the case of PG, is difficult to define in practice, but 
can be attributed to the presence or absence of applied stress. 

There are better data for growth of PG than there are for creep, 
although deformation under creep would be expected to exceed that of growth. 

In any case, creep of PG might cause important strains (and resultant 
stresses) in free standing chambers, especially at high temperatures. 
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7. Operating Stresses 

During motor firing, the inside surface of a free standing 
pyrolytic graphite thrust chamber will rise to within a few hundred degrees 
of the combustion temperature in areas of high heating rates such as the 
throat, while a large temperature gradient will exist across the wall, due to 
the low thermal conductivity in the c-direction. The operating stresses in 
the thrust chamber will consist of the total of four components: (l) pres¬ 
surization stresses, (2) residual stresses, ( 3 ) thermal stresses, and (4) 
creep and growth stresses. A complete structural analysis of the PG chamber 
requires analysis of the operating stresses for the complete history of 
the chamber. Thermal gradic its will change rapidly immediately after ignition, 
will reach steady state values, and then will change again following shutdown. 
Stresses may change even during steady state (temperature) conditions if the 
wall temperature is high enough to cause creep and growth. The stresses during 
a multiple firing will depend on the temperature history of the chamber and 
the length of elapsed time during previous firings. 

8 . Postrun Stresses 


If permanent deformation of the PG occurs due to growth, 
plastic behavior, or creep, the stresses in the chamber will differ from the 
as-deposited residual stresses. The residual stresses at room temperature, 
after one or more motor firing cycles, are defined as postrun stresses. 

C. Heat Transfer Problems 


Most of the important heat transfer problems in a free standing PG 
thrust chamber are caused by the very low thermal conductivity of the PG in 
the radial or c-direction. As a consequence, the heat sink effectiveness of 
ithe PG wall is small, and, if radiation cooling is the only cooling process 
in effect, the inner surface temperature will rise rapidly in regions of 
high heating rates such as the nozzle throat. Steady state temperatures at 
the inner surface will be reached within a few seconds, depending on com¬ 
bustion chamber conditions, and will be close to the combustion gas tempera¬ 
ture. 


Another heat transfer problem is the heating of the injector during 
and after motor operation. The high a-direction thermal conductivity will 
speed the conduction of heat along the chamber wall to the injector joint 
and to the exit nozzle. At the same time, the inside surface will act almost 
as a black body radiating to the injector face. During motor operation, the 
injector face will also be heated by convection and radiation from the com¬ 
bustion gas. This heat can be absorbed during firing by propellant cooling. 
After shutdown, the hot chamber walls will continue to radiate heat to the 
injector face. The inner wall temperature, (and hence the radiation flux) 
will decrease with time as the stored heat is conducted and radiated away. 

The principal problem caused by overheating the injector is that 
of propellant vaporization, causing ignition delay or pressure transients 
on restart. Injector heating after shutdown may also cause thermal stresses 
and structural problems in the injector, or may even cause melting of low 
temperature components or brazed joints in the injector. 
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D. Chemical Erosion 


•The combustion temperatures of the Fg/Hg and Fg/NgH^ blend will 
range between 5500* and 7500*F, depending on mixture ratio, combustion ef¬ 
ficiency, and chamber pressure. The inside surface of the PG chamber may 
approach the combustion temperature for long duration runs, raising the 
question of the rate of chemical erosion of the PG by the combustion gases. 
High erosion rates, coupled with the limitations on wall thickness, might 
severly restrict the run times that free standing PG chambers could endure 
without weakening the wall to the point of failure. An important point here 
is the fact that the rate of chemical erosion would be about the same for 
different size motors at the same chamber pressure. Larger chambers can 
use thicker walls, hence the allowable run times would be greater. The major 
combustion product of Fg/Hg and Fg/Hydrazine blend is hydrogen fluoride. 

Based on thermodynamic equilibrium calculations, it is predicted that HF will 
begin to react with graphite at temperatures above 4500*F. Other combustion 
products, such as hydrogen, nitrogen, and fluorine, will also react with 
graphite. However, the rate of chemical erosion is impossible to predict 
without experimental kinetic data, which are presently lacking. No analysis 
of the chemical erosion problem has been done in the current analysis and 
preliminary design phase. 

E. Analysis and Preliminary Design 

This Phase I report presents the results of a six-month analytical 
and design study which has been conducted to investigate the structural and 
heat transfer characteristics of free standing PG chambers. The analyses 
were performed for chamber pressures up to 300 psia. The following combina¬ 
tions of propellants and thrust levels were considered: 


1 . 

100 

-lb Thrust Level 


a. 

n 2 o 4 /o. 5 N 2 h 4 . 0.5 otmh 


b. 

Fg/BAlOlU* 


c . 

V H 2 

2. 

1000-lb Thrust Level 


a. 

NgO^/O^NgH^ - 0.5 UDMH 

3. 

5000-lb Thrust Level 


a. 

Fg/BA10l4* 


A pyrolytic graphite product evaluation was also made, consisting 
of a vendor survey of current production techniques, together with laboratory 
testing of 2-inch diameter PG tubes. 


* Hydrazine blend: 67 $ NgH^ + 24$ MMH + 9$ HgO 
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II. PRELIMINARY DESIGN STUDIES 


A. Scope of Studies 

Preliminary design studies were made to establish the size and 
shape of PG chambers for stress and heat transfer analyses. These design 
studies were made for those combinations of thrust and propellants which 
were to be used later in thrust chamber test firings, as follows: 

1. 100-pound Thrust (Sea Level ) 

N2O4/0.5 N2H4 - 0.5 UDMH 
F 2 /BA1014 

F2/GH2 

2. 1000-pound Thrust (Space) 

N2O4/O.5 N 2 H 4 - 0.5 UDMH 

3» 3000-pound Thrust (Space) 

F 2 /bA1014 


B. Propellant Performance 

The performance of the three propellant combinations which were 
studied is given in Table I. The design mixture ratios used for most of the 
analytical parameter studies were chosen as 2.0 for N2O4/O.5 N2H4 - 0.5 UDMH 
and F2/BAIOI4, and 12.0 for F2/H2* Most of the heat transfer and stress 1 
studies were done for thrust chamber sizes based on a C* efficiency of 95^, 
using shifting equilibrium. 

The variations in combustion temperature and specific impulse with 
changes in mixture ratio are given for N2O4/O.5 N2 H 4 - 0.5 UDMH in Figure 2, 
for F2/BAIOI4 in Figure 3, and for F2/H2 in Figure 4. 

The variation in combustion temperature with changes in C* 
efficiency is given for the three propellant combinations in Figure 5. 


C. Thrust Chamber Dimensions 


1. 100-pound Thrust Configuration 

Analytical parameter studies for the 100-pound thrust level were 
to be made for three different propellant combinations and for chamber 
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pressures of from 50 to 300 psia. Due to the large number of conditions to be 
studied, It was decided to choose a standard 100-pound thrust configuration 
and to study the effects of various propellant combinations and combustion 
conditions on that configuration. The parameters to be varied included the . 
vail thickness-to-radius ratio (t/r), which is a particularly significant 
parameter for pyrolytic graphite chambers, since both pressure stresses and 
residual stresses are related to t/r. Thus, much of the optimization of 
chamber dimensions is reflected in the effect of t/r. 

The 100-?ound thrust configuration is shown in Figure 6. This 
configuration had been previously manufactured and tested by several agencies 
and FG chambers were immediately available for further testing if desired. 

A comparison of the thrusts for the three propellam, combina¬ 
tions using the 100-pound thrust configuration is shown in Table II. It is 
interesting to note that the thrust would be almost identical at a given 
chamber pressure for all three propellant combinations. 

Wall thicknesses in the throat and combustion chamber of the 
100-pound thrust configuration are shown in Table III for various values of 
t/r. 


2. 1000-and 3000-pound Thrust Configurations 

The dimensions for the 1000-pound space thrust configuration 
using Np0^/0.5 NpH^ - 0.5 UDMH and for the 5000-pound space thrust configura¬ 
tion using F /BAlOlU are given in Table IV, using the nomenclature for dimen¬ 
sions given In Figure 7 . The dimensions of the 1000-and 5000-pound thrust 
chambers were varied with chamber pressure to keep the thrust constant. 

The chambers were sized using a radius of curvature (Rjj) at 
the throat of twice the throat diameter (D t ) and a radius of curvature (R) 
at the transition of the nozzle convergent section which was equal to the 
combustion chamber radius. The chamber lengths listed in Table IV are for 
a chamber L* of 40 inches. The nozzle lengths are for a 15 degree half¬ 
angle exit cone. 

The wall thicknesses for the 1000-pound and 5000-pound thrust 
configurations are shown in Tables V and VI, respectively. 
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III. PROPERTIES OF PYROLYTIC GRAPHITE 


A search of the literature was made to determine the most applicable 
thermal and mechanical properties of pyrolytic graphite in the temperature 
range of 0® to 7000*F. The properties of PG which were chosen to be used in 
the heat transfer and stress studies are called the "analysis standard" 
properties. In no case were any data found for temperatures above 5000 # F, 
but the analysis standard properties to be used for the parametric stress 
analyses are all extrapolated to 7500°F to facilitate computer operation during 
iterative calculations. 

PG will experience sublimation at temperatures near 6600°F. The 
flame temperature of F 2 /H 2 can be as high as 7500*F, but for the mixture ratios 
and C* efficiencies of practical importance, the wall temperature will not 
exceed 6500°F, and will be even lower if injector film cooling is used. 
Furthermore, the endothermic or exothermic effects of chemical reactions 
between the PG and combustion gases will influence the actual temperature of 
the inner surface of the thrust chamber. However, these effects cannot be 
calculated at this time since the kinetics of these reactions are not known. 

The analysis standard curves were established in two steps. First, 
all available data on a particular property was plotted on one graph. Then, 
a 'best" single curve was determined from this collection, taking into account 
when the data were measured. That is, more credence was given to recent data. 

No attempt was made to identify the differences in material properties from 
one vendor to another for several reasons. First, the material properties are 
undoubtedly dependent on microstructure. Second, a wide variety of micro¬ 
structures have been and are being made by each of the vendors so it is 
misleading to attempt to label some particular type of PG as being representative 
of one vendor. Finally, much of the past testing has indiscriminately mixed 
the results from different types of microstructure, and the effort required 
to review all of the' old data in detail was beyond the scope of this program. 

All that was attempted in drafting the analysis standard properties was to 
choose properties that probably were representative of much of the material 
currently being produced. 

Other than room temperature values, no data were available for the 
modulus of elasticity in the c-direction nor for the Poisson's ratios. There¬ 
fore, constant values of 0.95 for an< * ‘0»17 f° r v i2 were Poisson's 

ratios were assumed constant for all temperatures, the modulus in the c-direc- 
tion was assumed to vary in the same proportion as the modulus in the a-direc- 
tion, using the ratio of the room temperature values (E c /E a )Kptas a scaling 
factor. Although this arbitrary estimate of c-direction modulus seems to be 
highly questionable, it can be justified by the fact that the value of c- 
direction modulus does not heavily influence the stresses in free standing PG. 
Since only limited data were found for the compressive modulus, these data 
were included with the tensile data to give one modulus curve. Thermal 
expansion coefficients were determined as mean values from 0°F to the temper¬ 
ature of interest. The elastic constants at elevated temperatures were checked 
against certain restrictions derived from the strain energy and compressibility 
criteria given in Reference 5« These restrictions were not violated in any 
of the stress analyses. 
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The total normal emissivity in the c-direetion was taken to be 
constant at a value of 0.8. This was based on personal communication with 
R. J. Champetier, Aerospace Corporation. 

An extensive search was made for experiment ad data on the creep of 
FG in the a-direction. Several types of data were found (References 6 to 10). 
Only Kotlensky’s data (References 6 and 10) could be used for stress analysis, 
and these data were limited to creep at one stress level. The other refer¬ 
ences (References 7, 8, and 9) present either flexural creep data, which are 
not usable, or present exploratory data with temperatures varying during 
loading so that creep under constant loading and temperature could not be 
evaluated. A DDC literature search and contact with those organizations 
which have tested FG in the past shoved that no other creep data were avail¬ 
able. The results of the stress analyses showed that creep will have a large 
effect on the structural capability of free standing FG chambers for the high 
temperature propellants such as F2/BAIOI4. Therefore, more creep and growth 
data are needed. 

•The analysis standard curves are presented in Figures 8 through 13. 
The data scatter for two important properties, a-direction modulus and c- 
direction thermal conductivity, are shown in Figures 14 and 15. The data for 
Figures 14 and 15 were obtained from References 8 and 11 to 15. Again, it 
should be emphasized that the properties from each reference typically repre¬ 
sent materials from many sources. For example, the FG properties from Lock¬ 
heed (Reference 8) were actually measured by Aerojet General Corporation and 
Jet Propulsion Laboratory (JPL). Furthermore, the FG used in the tests was 
made by General Electric Co. (Schenectady), High Temperature Materials, Inc., 
General Electric Go., (Detroit), and Raytheon. One interesting comparison 
of the effect of microstructure is shown in Figure 14, which shows modulus 
of elasticity in the a-direction for substrate nucleated PG made by General 
Electric Co. (Detroit) and for continuously nucleated PG made by Raytheon Co. 

The data scatter for c-direction thermal conductivity (Figure 15) 
is greater than for any other property of FG. The principal cause is pro¬ 
bably due to the effects of delaminations and microcracks in the FG, as 
discussed in Reference 13. The trend is toward lower c-direction con¬ 
ductivity in FG which is free of microcracks and delaminations. 

The growth and creep data used in the stress analysis are plotted in 
Figure 16. The creep curve is taken from Reference 6 and it represents the 
uniaxial tensile creep of PG in the a-direction. The applied stress was 
12,000 psi and the temperature was 4712°F. The creep curve shown in Figure l6 
was one of several curves presented in Reference 6, all between 4200®F and 
5000*F. These curves showed a somewhat mixed trend for the 400 minutes 
duration of the tests, and they are also difficult to interpret at short times. 
The creep-versus-time curve lies above the growth-versus-time curve at 4712 # F, 
as expected. 
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IV. HEAT TRANSFER STUDIES 

Heat transfer calculations were made to determine the temperature 
gradients in free standing PG thrust chamber walls (required for thermal 
stress analysis) and the temperature rise of the injector after shutdown. 

A. Wall Temperature Gradients 

Transient heating calculations were made with an IB! 7040 computer 
program to evaluate the radial temperature gradients in the wall of free 
standing PG thrust chambers. The chambers were considered to be radiating 
to space from the outside surface, using an emissivity of 0.8 for the PG. 

The analysis standard thermal properties (See Section III) were 
used in all heat transfer calculations except where otherwise noted. The 
convective heating rates of the inside wall by the combustion gas were cal¬ 
culated from the Bartz equation as modified in Reference 16 to use the 
reference enthalpy method. Heating of the wall by radiation from the com¬ 
bustion gases was not considered. Two reasons Justifying this simplification 
were as follows: (l) Radiation heating is usually much less than convective 
heating unless the chamber dimensions are very large and (2) The inner PG 
wall will in many cases approach the combustion temperature, so that little 
change in steady state wall temperatures will result from an Increased 
heating rate. The transient heating rates would be somewhat underestimated 
by neglecting radiation, but in any case the convective heating rates cannot 
be predicted with any high degree of accuracy. Only radial heat conduction 
was assumed. 

All heat transfer and stress analyses presented in this report 
are for a C* efficiency of 95# and the design mixture ratio (See Table II) 
unless otherwise specified. 

1. 100-pound Thrust, N2O4/O.5 ®MH 


The transient temperature rise of the inside and outside sur¬ 
faces of the throat of the 100-pound thrust chamber with ^0)^/0.5 N^^-O^ 
UIMH are shown in Figures 17, 18, and 19 for chamber pressures of 50, 100, 
and 300 psia, respectively. Corresponding curves for the combustion chamber 
are shown in Figures 20, 21, and 22. The dimensions of the 100-pound thrust 
configuration are shown in Figure 6. The throat diameter is 0.75 in * and 
the chamber diameter is 2.0 ins. Various values of wall thickness (See 
Table III) were studied, since t/r (the thickness-to-radius ratio) is a sig¬ 
nificant parameter both for residual and pressurization stresses. It can be 
seen that the inner wall temperature is within a few hundred degrees of the 
combustion gas temperature, except for the thinnest walls and the lowest 
chamber pressure, where radiation cooling reduces the wall temperature some¬ 
what. 

Typical transient temperature gradients in the chamber wall 
can be obtained from the heating curves (Figure 23) for the combustion cham¬ 
ber at 100 psia and a wall thickness of 0.045 in. 
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a* Effect of C* Efficiency 

The effect of varying C* efficiency on the steady state 
Inner combustion chamber vail temperature is shown in Figure 24 for a chamber 
pressure of 100 ptla, an 0/F ratio of 2,0, and a vail thickness of 0.045 in. 

b. Effect of Mixture Ratio 


The effect of varying the mixture ratio on the steady state 
inner combustion chamber wall temperature is shown in Figure 25 . 

2. 100-pound Thrust, Fg/BA10l4 


The transient temperature rise of the inner and outer wall sur¬ 
faces of the throat and the chamber wall for the 100 -pound thrust chamber are 
shown in Figures 26 through 31* All of these calculations were for Fp/BA10l4 
at an 0/F ratio of 2.0 and a C# efficiency of 95#* 

Typical temperature gradients in the combustion chamber wall 
with F 2 /BA10l4 at 100 psia are shown in Figure 32. 

a. Effect of C* Efficiency 

The effect of varying the C* efficiency from 85 $ to 100 $ 
on the steady state inner combustion chamber wall temperatures is shown in 
Figure 33. 


b. Effect of Mixture Ratio 


The effect of varying the mixture ratio on the steady state 
combustion chamber temperature is shown in Figure 34. 

3. 100 -pound Thrust, F^Hg 


A comparison of the temperature rise of the throat of the 100- 
pound chamber using F 2 /H 2 at an 0/F ratio of 12 and F2/BA1014 at an o/F 
ratio of 2.0 is shown in Figure 35* A similar comparison for the combustion 
chamber i* shown in Figure 3 6 . It is seen that the difference between the 
wall temperatures for the two propellant combinations is very small. There¬ 
fore, all heat transfer and stress analyses for Fg/BAIO^ are also applicable 
to Fg/Hg. 

a. Effect of Mixture Ratio 


The effect of varying the mixture ratio of F 2 /H 2 from 10 
to 14 is shown in Figure 37* The wall temperatures vary by 250 # F at most 
over this range. 

4. Thermal Conductivity Effects in 100-pound Thrust Chambers 

The effectiveness of increased radial thermal conductivity in 
reducing the inner wall temperatures was studied by comparing the following 
materials: 
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1. Pyrolytic Graphite 

Thermal conductivity (c-direction) ~ 0.75 Btu/hr ft°F 
Specific gravity = 2«2 

2. Carb-i-tex 100 

Thermal conductivity (Against grain) = 3.0 Btu/hr ft°F 

Specific gravity = 1.38 

3. Carbon Fibre Composite 

Thermal conductivity (Against grain) * 10.0 Btu/hr ft # F 

Specific gravity = 1.38 

It vas assumed that the Carb-i-tex and composite chambers 
were formed so that the radial direction was against the grain. The tempera¬ 
ture variation of the thermal conductivity for pyrolytic graphite was in¬ 
cluded, but the conductivities for Carb-i-tex and the composite material 
were assumed to be constant. The specific heat for graphite (See Figure 10) 
was used for all three materials. 

The comparison was made for the throat of the 100-pound thrust 
configuration at 100 psia with a t/r of 0.12. The results for N2O4/O.5 
N2H4-O.5 UTMH (Figure 38) and for F2/BAIOI4 (Figure 39) show that a four-fold 
Increase in the radial thermal conductivity (from 0.75 to 3*0 Btu/hr ft # F) 
would reduce the steady state inner wall temperature by only 300 # F at most. 
There is a crossover in the transient heating curves because of the chang¬ 
ing importance of density and radial thermal conductivity as the wall tem¬ 
perature rises. 

5. 1000-pound Thrust, N2O4/5O-5O 


In contrast to the 100-pound thrust chamber analysis, which 
studied the same size chamber at different values of wall thickness and 
chamber pressure, the analyses of the 1000- and 5000-pound thrust configu¬ 
rations was done by changing the chamber size at various chamber pressures 
so as to keep the same thrust. The relationships of wall thickness and 
t/r for the 1000-pound thrust configurations is shown in Table V. 

The temperature rise of the 1000-pound thrust combustion 
chamber wall is shown in Figures hO , 4l, and 42 for chamber pressures of 
50, 100, and 300 psia, respectively. The temperature rise of the throat 
at 100 psia is shown in Figure 43. Heat transfer and stress analyses for 
the throat were made only at 100 psia, since additional analysis at 50 or 
300 psia was not considered warranted in view of the uncertainty in throat 
residual stresses. All of these curves are for N2O4/O.5 N2H4-O.5 U1MH 
at an 0/F ratio of 2.0 and a C* efficiency of 95#* 
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6 . 5000-pound Thrust, F 2 /BA10l4 


The temperature rise of the 5000 -pound thrust combustion cham¬ 
ber vail Is shown In Figures 44, 45, and 46 for chamber pressures of 50, 100, 
and 300 psla, respectively. The temperature rise of the throat at 100 psla 
Is shown In Figure 47. 

The heating rates for the 5000 -pound thrust chamber were based 
on using Fp/BA10l4 at a mixture ratio of 2.0 and a C* efficiency of 95$* The 
relationships of wall thickness and t/r for the 5000 -pound thrust chambers 
are shewn In Table VI. 

7* Controlled Wall Temperature Chamber 

Pyrolytic graphite chambers with a 4500°F temperature at the 
Inner surface were of Interest because thermochemical calculations show 
that 4500*F is the threshold for chemical erosion of pyrolytic graphite by 
HF, one of the principal combustion products of F 2 /BA10l4 or F 2 /H 2 , and 
because growth is negligible below 4500*F, according to the data in Figure 
16. Therefore, some stress and heat transfer analyses were made for such 
a chamber, referred to as the controlled wall temperature chamber. 

The steady state outside wall temperatures of this chamber 
for vcpious wall thicknesses are shown in Figure 48. An average c-direction 
thermal conductivity of 0.75 Btu/hr ft°F was used. The outer wall was 
radiation cooled with an emissivlty of 0 . 8 . 

B. Postrun Injector Temperatures 

Free standing PG chambers have a low thermal conductivity across 
the wall, and as a result the inner wall temperature is quite hot and cools 
down slowly, compared to most other types of thrust chambers. This raises 
the problem of possible overheating of the injector by radiation and con¬ 
duction from the PG walls after motor shutdown. This postrun injector 
heating is sometimes referred to as "soak back". 

1. Preliminary Soak Back Analysis 

A preliminary analysis of the postrun injector temperature 
rise (soak back) was made by studying the following simplified thermal models 
of free standing PG chambers: 

1. 100-pound thrust, Fg/BA10l4, P c = 100 psia 

2. 1000 -pound thrust, RgO^/O^ NgH^-O^ UEMH, P c = 100 psia 

3 . 5000-pound thrust, F 2 /BA1014, P c = 100 psia 

A sketch of the dimensions of those thermal models is shown in 
Figure 49. The injectors were assumed to be i/2-inch thick monel with an 
emissivlty of O. 25 . A check of these model injector weights with the actual 
weight of two Marquardt injectors (a 100-pound and a 1000-pound injector) 
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showed an adequate correlation, so that the model injector weights are con¬ 
sidered representative. 

The preliminary soak back calculations were made by assuming 
that the PG vails at the time of motor shutdown were at the combustion cham¬ 
ber steady state temperatures. Two cases were studied for each model: (l) no 
thermal resistance between chamber walls and injector and (2) infinite thermal 
resistance between chamber walls and injector. No heat loss from the nozzle 
throat or from the injector was considered in these preliminary analyses. 

The maximum injector temperatures reached during soak back are 
shown in Table VH. The temperatures of the chamber and injector after shut- 
dc -1 are shown for the 100-pound thrust chamber in Figure 50 for the case 
of no thermal resistance between the chamber and the injector. 

Similar data for the case of infinite thermal resistance between 
the 100-pound thrust chamber and the injector are shown in Figure 51* For 
the 100-pound thrust chamber, infinite (i.every large) resistance can 
reduce the injector temperature to an acceptable value of about 230°F. How¬ 
ever, as shown in Table VII, the soak back problem becomes more severe as 
the chamber size and wall thickness are increased. Therefore, a more detailed 
thermal analysis was made for the 5000-pound thrust chamber, as described 
below. 


2. Final Soak Back Analysis of 5000-pound Thrust Chambers 

The chamber contours for the final analysis of the soak back 
problem for the 5000-pound thrust chamber were more representative of an 
actual thrust chamber than the thermal model used in the preliminary soak 
back analysis. The chamber dimensions were those listed in Table VI, using 
the nomenclature of Figure 7* A summary of the analysis parameters is given 
in Table VIII. 


The analysis Included radiant reflections within the chamber 
and radiation losses through the throat and it was performed with the IBM 
7040 Thermal Analyzer Computer Program, using the radiosity method (Refer¬ 
ence 17). The emissivity of the inside surface of the PG chamber was taken 
as 0.80. Infinite resistance to conduction between the chamber walls and 
injector was assumed, since radiation from the inside wall surface to the 
injector is the principal mode of heat soak back for the 5000-pound thrust 
chamber. 


Heat soak back can be alleviated by increasing the injector 
weight, but it was desired that the effects of PG thrust chamber parameters 
(such as thickness) be studied separately. Therefore, the injector weight 
for the final 5000-pound thrust chamber soak back analysis was kept constant 
at 17.05 pounds. The initial PG vail temperature throughout the chamber and 
throat region was taken as the steady state temperature of the combustion 
chamber. The wall thickness throughout the chamber and throat region was 
taken as the wall thickness in the combustion chamber for the chosen con¬ 
traction ratio and a combustion chamber t/r ratio of 0.0^5« 
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Several assumptions were made for the initial PG wall tempera¬ 
ture to account for the possible effects of film cooling* 

a. Injector Eaissivity Effects 

The injector soak hack temperature could he reduced by 
lowering the emisrivity of the injector, since radiation from the PG walls to 
the injector is a large part of the total heat soak hack. The following 
values of injector emissivity were studied: 

1. Emissivity * 0.05 (Silver plated injector) 

2. Emissivity = 0.25 (Monel injector) 

3* Emissivity = 0.8 

The injector soak hack temperatures for these values of 
injector emissivity are shown in Figure 5~ for an initial inside PG wall 
temperature of 4500®F, which would require film cooling of the wall. Soak 
hack temperatures for an emissivity of 0.05 and an initial inside wall temper¬ 
ature of 5600°F are also shown in Figure 52. Decreasing the injector emissi¬ 
vity from 0.25 to 0.05 hy silver coating or polishing would reduce the injector 
temperature from about 700° to about 200 8 F. Available data indicate that 
silver does not react with fluorine up to about 700°F, and therefore it is 
possible that a silver coated injector might maintain a low emissivity even 
in combustion chamber conditions. However, examination of Monel injectors 
after test firings has shown them to be covered with deposits. The source 
and composition of these deposits, or whether they would occur in space opera¬ 
tion, are not known. Any deposits on the injector would increase its emis¬ 
sivity (perhaps to 0.8). With a high injector emissivity, the injector 
temperature would approach the inside wall temperature in a relatively short 
time, as shown in Figure 5 2 « Therefore, rapid cool down of the chamber wall 
is desirable. 


b. Wall Thickness Effects 


The rate of cool down of the chamber wall can be increased 
by decreasing the chamber wall thickness. A thinner wall cools down faster 
because of its lower heat capacity and reduced conduction resistance across 
the wall, as shown in Figure 53. 

c. Contraction Ratio Effect 


If the wall thickness is reduced, while keeping a constant 
chamber diameter, the t/r ratio of the wall will decrease, with possible loss 
of structural strength. However, if the chamber diameter (i.e., contraction 
ratio) is decreased in proportion to the decrease of wall thickness, the 
structural strength, as a first approximation, may be about the same, and the 
thinner wall should reduce the injector soak back temperature. In addition, 
the injector for smaller contraction ratios will have a smaller injector face 
area exposed to the radiation from the PG. 
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Injector soak back temperatures for the 5000-pound thrust 
chamber, for a constant t/r ratio of 0.045, are shewn for contraction ratios 
of 2, 4, and 6 in Figure 54. The injector soak back temperature could be 
reduced from about 700 # F to about 350°F by decreasing the contraction ratio 
from 4 to 2. The corresponding reduction in vail thickness would be from 
0.264 to 0.186 inch. 

d. Film Cooling Effects 


A lover inside vail temperature at shutdown would result 
in a lover injector soak back temperature. Therefore, the following cases of 
film cooling of the PG wall were analyzed for a 5000-pound thrust chamber, 
with a contraction ratio of 4.0, and injector emissivity of 0.25, and a com¬ 
bustion chamber t/r ratio of 0.045: 

1. Inside wall temperature at shutdown of 5600°F 
(No film cooling) 

2. Inside wall temperature at shutdown of 4500°F 

3. Inside wall temperature at shutdown of 2500°F within 
5 inches of the injector and 4500°F elsewhere. 

The effect of film cooling on injector soak back tempera¬ 
ture is shown in Figure 55> which shows a reduction from 1000°F, without film 
cooling, to 550°F for the case of a film cooled temperature of 2500°F to 
4500°F. 


e. Injector Cover 

Another method studied for reducing the injector soak back 
temperature utilized a pyrolytic graphite cover plate for the injector face, 
with appropriate openings for the injector holes. For the case of a 0.5-inch 
thick cover, a chamber t/r ratio of 0.045, an injector emissivity of 0.25, 
and an inner FG wall temperature at shutdown of 5600°F, the injector soak 
back temperature was 1100°F. In this case, the front face of the injector 
cover was assumed to be at a temperature of 5600°F at shutdown. The injector 
and the back face of the cases were assumed to be at ambient temperature at 
shutdown. The resulting high temperature of the injector was partly due to 
the amount of heat stored with' ~ the cover plate at shutdown. 

f. Conclusions 


It appears that a combination of smaller contraction ratios 
(or thinner PG walls), low injector emissivity, and film cooled PG walls will 
be required to limit the 5000-pound thrust injector soak back temperatures 
to acceptable values. 
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V. STRESS ANALYSES 

Stress analyses of pyrolytic graphite thrust chambers were made to 
investigate the capability of the free standing PG structure to withstand the 
thermal and pressurization stresses during steady state and cyclic operation. 
Most of the stress analysis was done using an IEM 70^0 computer program which 
was written for analysis of a cylindrical shell of pyrolytic graphite. The 
shape effects of the nozzle throat region were not included in the IB4 pro¬ 
gram. The variation of elastic properties with temperature was included in 
the analysis. Permanent deformation, or growth, of PG was also included. 
Plastic deformation and creep were not included in the analyses, except in 
a very approximate manner. The outside pressure was assumed to be zero, 
corresponding to operation in space. The method of analysis is described 
below: 


A. Method of Analysis 

1. Anisotropic Cylinder Analysis 

A program for the E3M 7040 computer was written for the gen¬ 
eral solution of the stresses in the center of an infinitely long cylindri¬ 
cal shell, using cylindrical coordinates r, 0, and z. As shown in Figure 1, 
the crystallographic c-direction of the pyrolytic graphite is parallel to the 
radial coordinate (r), and the crystallographic a-b plane of the pyrolytic 
graphite corresponds to a cylindrical shell containing the axial direction 
(z) and the circumferential direction (©). Internal pressure, temperature, 
and axial loads were assumed to be axisymmetric, resulting in plane axial 
strain vhereby all principal shear stresses and principal shear strains 
vanish. 


One of the main features of the program is its ability to handle 
temperature-dependent elastic and thermal properties. 

The technique used is to subdivide the cylinder into N concentric 
cylindrical elements, where the value of N may be adjusted to suit the varia¬ 
bility of the material properties (both elastic and thermal), the radial tem¬ 
perature gradients, and any other relevant conditions. Each of the N elements 
is then assumed to have homogeneous elastic and thermal properties, determined 
by the temperature in the center of the element, although the elastic proper¬ 
ties may vary from layer to layer. An expression for the stresses and dis¬ 
placements can be found for each layer and when the conditions of continuity 
at the internal boundaries as well as the axial load and conditions at the 
external boundaries are taken into account, the stresses and displacements 
can be found throughout the cylinder. 

Pyrolytic graphite is a transversely isotropic material with 
the axis of symmetry In the c-direction. Therefore, the generalized Hooke’s 
Law can be written as follows: 
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*zz - 8 12 a ee + 8 n a zz + 8 13 a rr + *0 (T - V + Y e M 

8 rr - fl 13 a 00 + 8 13 °zz + 8 33 a rr + «r (T " T D> + Y r M 

a ee “ c n *00 + c i2 f zz + c i3 c rr " x e “ V " T e ^ 

a zz B c 12 *00 + C 11 *zz + c 13 *rr “ X 0 “ V ‘ T 0 ^ 

a rr c C 13 *00 + °13 *zz + c 33 € rr " X r “ T D^ “ T r ^ 


Where 

CT 1J “ Stresses 
* i j = Strains 
s^ = Elastic compliances 
= Elastic stiffnesses 

cr^ = Coefficients of mean thermal expansion 

X t = Thermal stress coefficients 

= Growth strains 

= Growth stress factors 

T = Actual temperature 
T D = Datum temperature 

The system of symbols and subscripts is similar to that explained in Refer¬ 
ence 18 • 
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The pyrolytic graphite is assumed to he stress free at the 
datum temperature, which was assumed in the current study to be 0°F (room 
temperature, approximately). 

The stiffnesses can be written in terms of the compliances 

as follows: 
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The compliances can be defined in terms of the conventional 
elastic constants as follows: 
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Where 

■ Modulus of elasticity In the a-direction 

■ Modulus of elasticity In the c-direction 

Vjo “ Poisson's ratio between contraction In the b-airection and 
elongation In the loaded a-direction 

V 13 * ?oiaeon's ratio between contraction In the c-direction and 
^ elongation In the loaded a-direction 

The thermal stress coefficients (Aj) are related to the coeffi¬ 
cients of thermal expansion as follows: 

X 6 * (c ll + C 12> “8 + c 13 “r (17) 


2 C Of + C Qt 

13 0 33 r 


(18) 


In addition, a relationship exists between the growth strains and growth 
stress factors: 
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(20) 


For this problem, the only equilibrium equation not identically 
equal to zero is the following: 


d q rr °rr ~ q 98 
dr r 
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( 21 ) 
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The strain-displacement equations are as follows, where u = 
Radial displacement: 


rr 


‘60 


du 

3r 


( 22 ) 

(23) 


*zz = (constant) 


(24) 


By substituting these strain-displacement equations into the 
stress-strain equations and thence into the equilibrium equation, we obtain 
the following differential equation for the radial displacement (u): 


d 2 u 


1 du 
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(25) 


In the above equation, all elastic constants (c^, c^, c^, c^, t^) have 

been assumed to be independent of radius. That is, for the particular 
element we are concerned with, all material properties are assumed constant 
and equal to the values of these properties evaluated from conditions at the 
mid-point of the element. 

If we define A as follows: 


A 



(26) 


and define 0 as follows: 


(X r " V ( T ~ V * (°12 " C 3 * ( T r " V ^r_ dT (27) 

C 33 r C 33 dr 
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the solution to Equation (25) can be written as follows: 

u - (Cj^ + Vj) / + (Cg + V 2 ) r‘ A (28) 


where and Cg are constants of Integration and 


V 1 ■ S / r' A + 1 0 dr (29) 

V 2 - * 5J / / + 1 0 dr (30) 


This solution is not correct, however. if A = 1.0 or A = 2.0. It means that the 
solution, instead of containing r to various powers, has an exponential term 
In It. If, in the analysis, this constant should be equal to 1.0 or 2.0 for 
a particular element, an error message will be printed out giving the element 
containing the violation and for the sake of computation a value of A of 1.001 
or 2.001 will be assigned for that particular element. 

A linear temperature distribution of the form 

T - T D = a Q + a^r (31) 


is fitted for each element and Equation (28) for u (r) now becomes 


. . Ba + DC, + E (B + C) a. r c 

jA + c r -A + - 2—1 -r + --L-L_ 

2 1 - A 2 4 - A 2 


UNCLASSIFIED 

-24- 



AFRPL-TR-66-95 


UNCLASSIFIED 


Report 6106 


Where 


B = 


C = 


< X r - V 

C 33 

^ c 12 ~ c 13 ^ 
°33 
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(36) 


The equations for the stresses can now easily be found as follows: 
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The constant C, delends on the axial restraint. If the ends of 
the cylinder are assumed fixed, C 3 ■ 0. For a finite axial strain, C 3 will 
not be zero and must be determined from the boundary condition 



Where: 

F s Axial load 

zz 

* Inside radius of cylinder 
r = Outside radius of cylinder 
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Since tne expression for o will contain different numerical coefficients 
for each element, this integration must be broken into N parts, as follows: 


/ r dr = 


IS o,n 

J ! °z 2 ,„ rdr 

” = 1 r i,n 


Therefore, from Equations (38), (40), and (4l), 
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we get 
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Since there are two unknown constants for each of the N elements 
(Cl^ n and C2,n) P lus the unknown axial constant C3, we need a total of (2N + 1) 
equations to solve for these (2N + l) unknowns. Two equations can he obtained 
for each boundary between any two elements, expressing the fact that the radial 
stress and radial displacement are continuous functions, as follows: 
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This gives 2(N-l) equations, since there are (N-l) common boundaries. In 
addition, at the inner boundary a ■ - P , that is 

XT O 


a rr ll = - Pq (internal pressure) (53) 

r i,i 


and at the outer *• inclary = 0, that is 



(For operation in space) 


(54) 


These two equations and the axial load equation, Equation (50), when added 
to the previous 2(N-l) equrtions give a total of 2N + 1 equations. The 
computer program solves for the (2N+1) unknowns and then substitutes them 
into Equations (37)# (38)# (39)# and (32) to give the complete stress and 
radial displacement distribution. The stresses calculated in this manner 
have a saw-tooth distribution across the wall, which can be smoothed out 
by increasing the number of cylindrical elements. A maximum of thirty 
elements are accommodated in the current computer program. 

Room temperature is a convenient datum temperature because 
thermal expansion data are usually expressed as a mean coefficient, from 
room temperature to the actual temperature. If the datum temperature is 
not 0°F, a mean coefficient of thermal expansion must be used which is 
calculated es follows: 


or 


i 



t d 


Where 

cr_ = Mean coefficient of thermal expansion from room 

i R temperature T_ to actual temperature T 
R 

or = Mean coefficient of thermal expansion from room 
T D temperature T R to datum temperature 


If T~ = Tr, this equation reduces to or m 

DR, 
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Room temperature is also a convenient datum temperature because 
all experimentally determined residual stress data are collected at room 
temperature, so that subsequent superposition of the residual stress effects 
is simplified. 

2. Superposition of Residual Stress Effects 

The theory of elasticity, as used in this analysis, is based on 
the assumption that the material is stress-free in the initial state, before 
the application of temperature changes or mechanical loads. If the material 
actually contains initial stresses, the stress distribution after application 
of mechanical loads and temperature changes can be obtained, in the case of 
perfectly elastic materials, by superposition of the initial stresses on the 
incremental stresses (defined as the stresses caused by loading and tempera¬ 
ture changes). In fact, initial stresses generally do occur in pyrolytic 
graphite, and are referred to as residual stresses. 

If, as in the present analysis, the variation of elastic modulus 
with temperature is considered, the method of superposition cannot be used 
to get the total stresses, since the residual stresses are related to a room 
temperature value of elastic modulus, whereas the incremental stresses typi¬ 
cally are related to such elevated temperatures that the elastic modulus 
has changed. 


The procedure used for including the effects of residual stresses 
in this study was to superimpose, on the incremental stresses calculated by 
the theory of elasticity, a residual stress contribution calculated from the 
initial (residual) strains and the elevated temperature elastic properties 
at the point in question. 

Since the residual stresses as found experimentally in as- 
deposited PG cannot be predicted with any success by theoretical means, it 
was necessary to use experimental residual stress data. Unfortunately, even 
the experimental residual stress data are not available in sufficient quantity 
to permit accurate extrapolation to other thickness and shapes of PG than 
those from which the residual stress data are collected. This problem will 
be discussed in more detail in Section V-B. 

For superposition of the residual stress contribution, it was 
assumed that only the residual strains are constant and that the residual 
stress contribution would vary as the material properties varied. This is 
only an approximation since it does not consider the redistribution of re¬ 
sidual strains for the case of nonuniform elastic properties associated with 
radial thermal gradients. A further assumption was made that the radial 
residual stress is zero, which is justifiable since the magnitude of the 
radial residual stresses within the wall is found experimentally to be low. 

In general, the data to be used were in the form of residual axial and 
circumferential stresses on the inside and outside surfaces of cylinders 
where the radial stresses are in fact zero. A linear distribution of cir¬ 
cumferential and axial residual stresses across the wall was assumed. The 
residual strains were then calculated from the following stress-strain 
equations: 
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By writing « in terms of and t , and the substituting into Equations 
(58) and ( 59 )>the following Equations can be derived: 
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Values of the residual strains €*q and « z at the inner and 
outer surfaces were calculated from experimental residual stress data and 
input to the computer. Linear distributions of these residual strains aross 
the wall were assumed, and at any point in the wall the stress component 
caused by the residual stresses was calculated from Equations (60) and (6l), 
using the temperature-dependent elastic properties at the point in question. 
These residual stresses were then superimposed on the incremental stresses 
to get the total stresses. 

3« Growth Analysis 

The effect of growth of PG was included in the computer stress 
analysis by using values of circumferential growth strain (X 0 ) which were 
estimated for each element from the data of Reference 3, based on the tempera¬ 
ture and elapsed time. The effects of temperature variations within each 
cylindrical element were approximated in a manner analogous to strain harden¬ 
ing laws. 


The c-direction growth was assumed to be twice the magnitude, 
and opposite in sign, to the a-direction growth shown in Figure 16. The 
growth strains in both the axial and circumferential directions were the same. 
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4. Creep Analysis 

Creep effects were included in the computer analysis by using 
a value of growth strain (Xq) which included a contribution caused by creep. 
The amount of creep was estimated from the probable level of stress in each 
element. Many simplifications were necessary. For example, creep was as¬ 
sumed to be directly proportional to stress level, using linear interpolation 
of the data in Figure 1 6 , between the zero stress curve (growth) and the 
12,000 psi curve. In actuality, the creep of most materials is roughly de¬ 
scribed by the following equation: 


e 


c 


= A o 


n 


(62) 


Where 

e = Creep rate 
c 

a = Stress 
A = Constant 

n = Constant (Almost always greater than l) 


The values of n for pyrolytic graphite cannot be estimated from 
available data. It was also assumed that creep rates under compression in 
the a-direction were of the same magnitude, but opposite in sign, as the creep 
rates for tension. 

The creep at temperatures above 4700°F was estimated using scaling 
factors obtained from Figure l6. The computer program, desi'med only for 
growth analysis, could not distinguish between axial and circumferential creep 
due to the different stress levels in these two directions. Therefore, the 
estimated circumferential creep was also used for axial creep. 

A computer program which makes a more accurate analysis of creep 
is feasible, but was not undertaken because of the scarcity of creep data. 

5« Plastic Deformation 


The current analysis can consider the temperature dependency of 
the elastic modulus, but it cannot consider plastic deformation. The effects 
of plastic deformation will not be important until the stress level and temper¬ 
ature are quite high. For example, at 4500°F, the plastic a-direction strain 
at a stress of 20,000 psi will be only 0.1$, compared to an elastic strain of 
0.8$. Plastic strain will be more important at higher temperatures, but ade¬ 
quate data were not available for analysis. 
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B. Experimental Residual Stress Data 

Experimental residual stress data were collected In an attempt to 
permit prediction of residual stresses in PG thrust chambers of various sizes 
and vail thicknesses. Most of the data available were for axial and circum¬ 
ferential residual stresses at the inside and outside surfaces of cones and 
cylinders. More detailed data, including radial stresses and residual stresses 
within the wall, are reported for a few cases in References 1, 2, 4, and 19 . 
The results of these studies show that residual stresses in PG are a function 
of many parameters such as the deposition conditions and thickness. At the 
present time, although some interesting hypotheses have been suggested for 
explanation of as-deposited residual stresses, the subject still remains 
largely undefined. 

In the course of this work, the residual stresses on cylinders and 
thrust chambers were measured. The limited scope of this work was sufficient 
to show only that the extrapolation of cylindrical residual stress data to 
other thicknesses or shapes is not possible at the present time. 

1. Cylinders and Cone B 

Experimental values of axial and circumferential residual 
stresses on the inside and outside of free standing PG cylinders and cones 
are shown in Figure 56 . The open symbols are old data and the closed symbols 
are new data obtained during this program. The analysis standard curves for 
residual stresses to be used in the stress analysis are also shown in Figure 
56 . 


Of necessity, the analysis standard curves were used for thrust 
chamber sizes and wall thickness/radius ratios far beyond the range of 
available data. Some of the new data obtained in this program were cor¬ 
rected by the ratio of the modulus of elasticity of the materials (boron- 
doped and continuously nucleated) and that of average PG (assumed to be 
4.5 x 10° pel). The as-deposited extrapolations are in fairly good agree¬ 
ment with the linear extrapolations up to a t/r ratio of about O.O 85 . No 
experimental data are available from tubes at t/r ratios above O.O 85 . 

2. Residual Stresses in Thrust Chambers 


The residual stresses in a cylindrical or conical combustion 
chamber would be expected to be similar to the residual stresses measured 
in cones and cylinders. However, the residual stresses in the vicinity of 
the throat of a free standing PG thrust chamber are complicated by the fact 
that two curvatures, the throat radius and the radius of curvature, have an 
influence on the residual stresses. The axial stresses are particularly 
affected by this double curvature. The curvature at the beginning of the 
contraction region will also cause additional axial stresses. The method of 
analysis used in calculating these stresses was that reported in Reference 
20, wherein good correlation was obtained between one data point of experi¬ 
mental axial residual stress and an analysis model. This model for axial 
stresses at the throat consisted of superposition of the following: 
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1. As-deposited axial stresses, based on the analysis standard 
curves for tubes (Figure 56 ) and the cross sectional di¬ 
mensions at the thrust chamber section of interest. 

2. Circumferential residual stresses in a hypothetical cylinder 
with a radius equal to the radius of axial curvature at the 
thrust chamber section of interest. In the throat region, 
this hypothetical cylinder, assumed to be stress-free at 
4000°F, would have an inside surface coinciding with the 
outside surface of the thrust chamber, as shown in Figure 

57• At the transition (Section A) to the convergent section, 
the hypothetical cylinder would have an inside surface coin¬ 
ciding with the inside surface of the thrust chamber. 

The above model was tested by measuring the residual stresses in 
two FG chambers of the configuration shown in Figure 58 . 

A comparison of the measured residual stresses in the PG chambers 
and the predicted residual stresses is shown in Table IX. The residual stress¬ 
es were measured by strain gage isolation. Predictions were made both by 

use of the analytical model described above for consideration of the nozzle 
shape, and also on the basis of the analysis standard extrapolations for 
cylinders. Examination of Table IX leads to several conclusions, as follows: 

1. The circumferential residual stresses in the throat were 
only about 40$ as great as those predicted by the analysis 
standard curves for cylinders. 

2. The axial residual stresses on the outside of the throat 
were the highest stresses measured, and they were far greater 
than those predicted by either the cylindrical or the nozzle 
shape models. 

3. The residual stresses in Section A were more closely pre¬ 
dicted by the cylindrical model (ignoring the effects of 
axial curvature) than by the analysis model for axial curva¬ 
ture. This is probably due to the fact that the amount of 
curvature at Section A was actually rather small. 

The results of this comparison led to the conclusion that the 
reliability of axial residual stress predictions at the throat of PG thrust 
chambers is so low that analysis of these stresses was not warranted. There¬ 
fore, no operating or postrun social stresses at the throat are presented in 
this report. A second conclusion was that the analysis of stresses in chambers 
with large t/r ratios (say 0.14) is not very accurate, using the analysis 
standards for as-deposited residual stresses. It is therefore obvious that 
much more data are required to permit evaluation of the effects of residual 
stress on the total stress in free standing pyrolytic graphite thrust chambers. 
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C. Proof Stresses 


Proof stress Is defined as the pressurization stress superimposed 
on the residual stress. Room temperature is assumed. The circumferential 
proof stress in a cylinder (which is assumed to represent any cross section 
of the thrust chamber) is shown, for the inside surface, in Figure 59. It 
is shown that the optimum t/r ratio, to minimize the proof stress, is at very 
low t/r values (about 0.02) for an internal pressure of 100 psla, while shift¬ 
ing to higher values for higher pressures. 

The circumferential proof stress on the outside surface is shown in 
Figure 60, which shows steadily decreasing stress with increasing t/r ratio 
values. The stresses at the inside and outside surfaces are used in most of 
the stress analysis since the stresses are almost always maximum at these 
surfaces. 


It should be noted that the residual stresses and proof stresses 
discussed above will apply, in view of our assumptions, to all sizes of thrust 
chambers, since both the pressurization stresses and the as-deposited residual 
stresses are functions of the t/r ratio. The analysis of 100-pound thrust 
chambers was done using the basic C100-1 dimensions, that is, a of 0.75 
inch and a chamber diameter of 2.0 inches. Various wall thicknesses were 
used. The analyses of the 1000-pound and 5000-pound thrust chambers were 
made for various size chambers, keeping thrust constant at chamber pressures 
of 50, 100, and 300 psia. Unless otherwise specified, all analyses were for 
a C* efficiency of 95^, the design mixture ratio (Table I), and a chamber 
pressure of 100 psla. 

D. 100-pound Thrust Configuration, 1^0^/0.5l^H 1+ -0.5 UDMH 


Steady state operating stresses in the combustion chamber of the 
100-lb configuration using NgOy O. 5 N 2 H 4 -O .5 UDMH are shown in Figure 6l as 
a function of t/r ratio. A similar curve for the throat is shown in Figure 
62 . The wall temperatures of the 100-pound thrust chamber with IfeOl*/ 0.5Nk H 4 
-0.5 UDMH were not quite hot enough to cause growth, which, according to 
Reference 3, is zero for 10 minutes at 4700 # F. 

Steady state operating stresses at the throat of the 100-pound 
thrust chamber with N 2 O 4 / O. 5 N 2 H 4 -0.5 UDMH are plotted versus chamber pres¬ 
sure in Figure 63 . 

Comparison of the proof stresses (Figures 59 and 60) with the oper¬ 
ating stresses (Figures 6l, 62, and 63 ) shows that the proof stresses are 
always higher than the operating stresses. In either case, the maximum stress 
is about 8000 psi tension, which is well within the capabilities of FG. 

E. 100-pound Thrust Configuration, Fg/BA10l4 


1. C» Efficiency Effect 

The variation of steady state operating stresses in the • 
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combustion chamber with variations in C* efficiency was calculated for the 
100-pound thrust chamber for F 2 /BAIOI 4 , using a wall thickness/radius ratio 
of 0.045* The axial and circumferential stresses at the inner and outer 
surfaces are plotted versus C* efficiency in Figure 64 and it is shown that 
the effect is rather small, for C* efficiencies between 90 and 100$. Most 
of the stress analysis was done using a C* efficiency of 95$* The maximum 
wall temperature for this comparison was 4520°F for a C* efficiency of 100$* 
Therefore, creep and growth were not included. 

2. 0/f Ratio Effect 

The variation of steady state operating stresses in the 
F 2 /BA10l4 combustion chamber with variations in o/F ratio were calculated 
for a 100-pound thrust chamber with a t/r ratio of 0.045. The effect of 
o/F ratio on combustion chamber stresses is very small, as shown in Figure 
65 . Most of the stress analysis was done using an o/F ratio of 2.0. The 
maximum wall temperature for this comparison was 4510 # F for an o/F ratio of 
2.5. Therefore, growth and creep were not included. 

3. Operating Stresses 

The steady state operating stresses in the combustion chamber 
of the 100-pound thrust F 2 /BA10l4 chamber were calculated for wall thickness/ 
radius ratios of 0.03, 0.045, 0.06, 0.08, and 0.10. The steady state oper¬ 
ating stresses in the combustion chamber are plotted versus t/r ratio in 
Figure 66 . The maximum tensile stresses on the outside surface, (both axial 
and circumferential) decrease with increasing wall thickness, and did not 
exceed 8500 psi at any t/r ratio. 

The maximum temperature of the combustion chamber at 100 psia 
does not exceed 4600°F, so that growth will not affect the stresses. Some¬ 
what higher temperatures will occur at the throat and growth must then be 
considered. 


Transient stresses were calculated in the wall of the com¬ 
bustion chamber with a t/r ratio of 0.045 and a pressure of 100 psia, and 
they are shown in Figures 67 through 70, which show the axial and circum¬ 
ferential stresses at both the inside and outside surfaces. The stresses 
at the surfaces are nearly always greater than those within the wall. Exami¬ 
nation of these curves shows that the stresses soon after motor ignition, 
when the radial temperature gradients are highest, are less severe than 
either the proof stresses or the steady state operating stresses. Creep and 
growth were not included in these calculations, since the wall temperature 
did not exceed 4300°F. 

Steady state operating stresses at the throat of the 100-pound 
thrust chamber with F 2 /BA10l4 were calculated assuming that growth had oc¬ 
curred over a firing time of 600 seconds. The stresses at the end of this 
time are shown in Figure 71 as a function of t/r ratio. 

The importance of residual stresses in calculating operating 
stresses is illustrated in Figures 72, 73> 74, and 75, in which the steady 
state operating stresses for the 100-pound thrust chamber with F 2 /BA10l4 
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are compared vlth the incremental stresses* The incremental stresses are 
those stresses due to pressure and temperature and they represent the operat¬ 
ing stresses for a thrust chamber which is stress-free at room temperature, 
i*e*, no residual stress* As- discussed before, the total stresses are ob¬ 
tained by superimposing the incremental stresses vlth the residual stresses 
corrected to actual wall temperatures* This superposition procedure is an 
approximation and, in effect, superposes strains rather than stresses, which 
is preferable if the material is not perfectly elastic. It can be seen from 
Figures 72 through 75 that the residual stress contribution to the total 
stress is very important* 

An analysis of the effect of growth during 10 minutes of firing 
was made in the throat section of the 100-pound thrust chamber for F 2 /ba 1014. 
Creep was not considered. Plots of circumferential stress versus firing time 
for t/r ratios of 0.08, 0*12, and 0.l6 are shown in Figures 75, 76, and 77 . 
The maximum tensile stress occurring anywhere within the wall is shown, and 
it usually is the surface stress. The operating stresses are generally with¬ 
in the range of 10,000 psi tension to 12,000 pel compression for firing t?mes 
up to 600 seconds* The stress levels were outside the above limits for short 
times (less than 1 second) for thick walls with t/r ratios of 0.12 and 0 .l 6 . 
This was due to the analysis standard residual stresses being so high at these 
large values of t/r, thereby being the dominant stress at such short times. 
After steady state temperatures are reached, the changes in operating stresses 
are due to growth. 

4. postrun Stresses 


Circumferential postrun stresses in the throat of the 100 -pound 
thrust Fp/lJAlOlU chamber for t/r ratios of 0.08, 0.12, and 0.l6 are presented 
in Figures 78 , 79, and 80. Growth was included in this analysis, but creep 
was not. Therefore, these plots of postrun stresses can be compared with the 
operating stresses in Figures 75, 76 , and 77 . 

In general, growth of the throat of the 100-pound thrust F 2 /BA10l4 
chamber does not cause operating or postrun stresses much greater than the 
original residual stresses, for firing times up to 10 minutes. However, it 
must be remembered that creep was not included in the above analysis. 

F. 100-pound Thrust Configuration - Comparison of Propellant Combinations 

The steady state operating stresses for the 100-pound chamber were 
compared for the two different propellant combinations and are shown in Figure 
8 l. Only the outside surface stresses are shown, since these are the maximum 
stresses. It can be seen that the stresses resulting from using the NgO^/ 
O. 5 N 2 H 4 -0.5 UDMH propellants are consistently lower than the stresses re¬ 
sulting from using F 2 /bA 1014. This is also true for the inside surface stress¬ 
es and is due to the'lower temperatures reached with NgO^/0.5NgH^-0.5 UIMH. 

G. 1000-pound Thrust Configuration, 1^0^/0.51^^-0.5 UDMH 


The steady state operating stresses in the combustion chamber of 
the 1000-pound thrust chamber are shown in Figures 82, 83 , and 84 for chamber 
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pressures of 50, 100, and 300 psia. The propellants were 1^0^/0.5I^Hi | -0*5 UDMH 
at a C* efficiency of 95 $ and a mixture ratio of 2,0, The stresses decrease 
for increasing t/r ratio tn a similar manner to those shown for the 100-pound 
thrust chamber in Figure 6l. The magnitude of the operating stresses is 
generally slightly higher than the 100-pound thrust chamber operating stresses, 
but they are still well within the stress limits of PG. 

H. 5000-pound Thrust Configuration, Fg/BA.1014 


1. Operating Stresses 

Operating stresses in the combustion chamber of the 5000-pound 
thrust chamber after 600 seconds of firing are shown in Figure 85 for a 
chamber pressure of 100 psia. The propellants were F 2 /BA10l4 at a C* effi¬ 
ciency of 95$ and a mixture ratio of 2.0. The stresses are fairly constant 
over the short t/r ratio range considered. Growth was considered for this 
t/r parameter curve but creep was not. 

An analysis was made of the effect of creep of PG on the oper¬ 
ating and postrun stresses of the 5000-pound thrust chamber at a chamber 
pressure of 100 psia and a t/r ratio of 0.045» 

The steady state inside surface temperature of the combustion 
chamber, using F 2 /BA10l4, at an o/F of 2.0 and 95$ C* efficiency, was 56ll°F. 

The inclusion of creep effects had a similar effect on both 
the circumferential and axial stresses, and, therefore, only the circum¬ 
ferential stresses will be discussed. Figure 86 shows the circumferential 
operating stresses versus firing time. 

The solid curve presents the stress analysis results without 
any creep considered. The dashed curve gives the results with creep included. 
Growth of PG was included in both analyses. By comparing the solid and dashed 
curves in Figure 86, it can be seen that creep changes the operating stresses 
by a very large amount. The stress on the outside surface was reduced by 
creep from 14,500 psi to 7500 psi after 600 seconds. The stress on the out¬ 
side surface was changed by creep from - 12,000 psi to + 7,000 psi after 
600 seconds. 


Figure 86 also shows the maximum tensile stress (circumfer¬ 
ential) occurring within the PG wall. As can be seen, this always occurs 
at the outside surface in the absence of creep, but sometimes occurs within 
the wall if creep is considered. 

2. Postrun Stresses 


The postrun stresses in the 5000-pound combustion chamber after 
various lengths of firing time are shown in Figure 87 . The postrun stresses 
are calculated assuming that the chamber has cooled down to room temperature. 
The internal stresses are now different in magnitude from the original re¬ 
sidual stresses because permanent deformation of the PG has occurred due to 
growth and creep. 
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The maximum postrun stresses on the inside surface are in¬ 
creased by the action of creep from - 32,000 psi (without creep) to + 1*0,000 
psi (with creep) after 600 seconds. 

The radial temperature drop across the 5000-pound thrust chamber 
wall was such that all creep and growth occurred within 0.100 inch of the 
inner surface. 


It is concluded that creep of PG has a very large effect on the 
stresses of chambers if the wall temperature is at least 4700°F. 

It appears that creep may improve the durability of a chamber 
for a single continuous firing, but would be damaging to the strength of the 
chamber for cyclic operation. The critical condition for cyclic operation 
would be motor ignition after cool down from a previous firing. It is also 
possible vhat a chamber might fail during cool down, due just to the internal 
postrun stresses exceeding the tensile strength of the PG. In any case, more 
creep data are required before precise conclusions can be drawn about stress¬ 
es in free standing PG chambers. 

I. Controlled Temperature Chamber 

Stress analysis was made of free standing PG chambers which were 
assumed to have an inner surface temperature of 4500°F. A linear tempera¬ 
ture gradient across the wall was assumed, with outer wall temperatures 
(for radiation cooling) as presented in Figure 48. 

The operating stresses on the inner and outer surfaces are pre¬ 
sented in Figure 88 for a 100-pound thrust chamber at 100 psia, and in 
Figure 89 for a 5000-pound thrust chamber, also at 100 psia. The trend for 
both sizes of chambers is to lower operating stresses at large values of 
t/r ratio. For example, at a t/r ratio of 0.10, all operating stresses 
would be less than 8000 psi, which is well within the strength capabilities 
of PG. 


Therefore, if the inner surface could be kept below 4500°F by film 
cooling, the PG chambers should be structurally sound. 

However, the calculated operating stresses depend heavily on the 
as-deposited residual stresses. In general, because of the anisotropic 
nature of PG, high residual stresses are reversed into low operating stress¬ 
es as the charrber temperature rises. Unfortunately, the available data on 
residual stresses at t/r ratios greater than 0.08 are not sufficient to give 
a firm foundation for stress analysis of thick walled chambers. This prob¬ 
lem is especially acute in the case of 5000-pound thrust chambers which will 
probably have thicker walls than ever before produced in cylindrical shapes. 

J. Expansion Nozzle 

Some analyses were made on the expansion nozzle of the 100-pound 
thrust altitude chamber at sections corresponding to expansion ratios of 
20:1 and 40:1. The operating stresses are due almost entirely to thermal 
and residual stresses. At the 20:1 section, the circumferential stress due 
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to pressure is only 21 psi and the axial stress due to axial load is merely 
7 psi. At the 40:1 section, there is no axial load and the circumferential 
stress due to pressure is only 13 psi. The total operating stresses for both 
sections are similar, because, while the thermal stresses are lower at the 
exit, the residual stresses are higher. At the inside surface, the circum¬ 
ferential stress is about 4300 psi compression, and the axial stress is about 
5900 psi compression. On the outer surface of the nozzle, the circumferential 
stress is about 3900 psi tension and the axial stress is about 6400 psi 
tension. It appears that the expansion nozzle of a free standing PG chamber 
will not be a problem area. 

K. Shear Stress 

The maximum shear stress (t) acting on the a-c plane in the a-direc- 
tion was calculated from the principal stresses ( a^ and a z ) using the 
formula 


T 


2 


o 

z 


No data were found regarding the shear strength on this plane. In all of the 
parameter studies, the operating shear stress never exceeded 3000 p6i except 
for several cases of postman stress. Sometimes PG tubes fail during pro¬ 
duction in a spiral failure, which might indicate shear failure. 

L. Vibration Stress Analysis 

A stress analysis was made of the 100-pound thrust chamber when 
exposed to vibration. The assumed power spectral density for the vibration 
loading is as follows: 

10 to 90 cps 


90 to 250 cps 
250 to 2000 cps 


0.055 g /cps at 10 cps, increases at 
3 decibles (db) per octave to 0.5 g^/cps 
at 90 cps 

2 

Constant at 0.5 g /cps 
2 

0.5 g /cps at 250 cps decreases at 
3 db/octave to 0.06 g^/cps at 2000 cps 


A preliminary analysis was made to determine the order of magni¬ 
tude of the stresses to be expected in the 100-pound thrust altitude chamber. 
The configuration analyzed was a 100-pound thrust chamber, except that a 
40:1 expansion nozzle was added. The critical mode of vibration was deter¬ 
mined to be the transverse mode, which causes maximum stresses due to bending 
moments from the transverse loading. 


UNCLASSIFIED 

-4l- 



1 


UNCLASSIFIED 


AFRPL-TR-66-95 


Report 6106 


To determine the maximum g level expected, the natural frequency 
and amplification factor of the chamber setup were v'-jtimated. The choice 
of natural frequency is unimportant because of the shape of the power spectral 
density distribution, and the amplification factor was taken to be 10, a 
typical value for many materials. The exact value of the amplification factor 
is not too important because it is raised to the l/2 power in the equation for 
the g level, as follows: 


■V? 


2 

(J-) 


cp6 r f f n< AF > 
n 


m 


Where 


2 

f n 
A F 


Power spectral density distribution value at f^ 

Natural frequency 
Amplification factor 


The 1 sigma g level is 44.3 g and the 3 sigma level is 133 g. This 
3 sigma level is the usual maximum level used in calculations, since only 1.2 
percent of all cycles are at or above this level. 

The vibration stresses in the PG chamber at the throat section and 
the injector flange were analyzed. The maximum vibration stresses (at the 
3 sigma level) are as follows: 


Location 

Shear 

(psi) 

Bending 

(psi) 

Injector flange 

150 

2,100 

Throat 

275 

5,000 


The vibration stresses are greater at the throat than at the in¬ 
jector flange, since the section properties of the throat are less, although 
the transverse bending moment is less. Although the bending stress of 5*000 
pel is appreciable, it is within the strength capabilities of good quality PG. 
The axial residual stresses at the throat would have to be added to the bend¬ 
ing stress due to vibration to get the total stress. 
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M. Tensile Strength 

Up until now, nothing has been said about the tensile strength of 
PG. Stress levels have been discussed without regard to possible failure 
of the chamber due to exceeding the tensile strength values. Curves of 
tensile strength of PG in the a-direction versus temperature are shown in 
Figure 90* It is seen that up to about 3000°F the strength is fairly constant 
and above 3000°F the tensile stiength increases very rapidly. This variation 
must be taken into account when checking operating stress Revels for failure., 
A stress level which will cause failure at room temperature may be well with¬ 
in the upper limit of tensile strength if the stress is reached when the 
chamber is hot. Considerations of biaxial loading or the effects of creep 
on the strength of PG cannot be evaluated at the present time. 

Work is being done to increase the room temperature tensile strength 
of PG. It is assumed that the same general behavior of increasing tensile 
strength above 3000°F will be maintained. 

In all cases examined in the parameter studies, the operating stres¬ 
ses were at no time above the average tensile strength corresponding to the 
temperature of the section where the stress is observed. 

For the postrun stresses, the allowable tensile strength was ex¬ 
ceeded in several cases. For high t/r ratios (greater than 0.12), the re¬ 
sidual stress exceeded the average tensile strength. If these predictions 
of residual stress were correct, as well as the tensile strength data, it 
would mean that tubes with t/r ratios greater than 0.12 cannot even be fab¬ 
ricated. This has been confirmed in production of tubes. However, thrust 
chamber throats have been fabricated with t/r ratios as high as 0.l6, for 
reasons not yet understood. 

For practically all cases where growth and creep were considered, 
the postrun stresses for firing times in excess of a few minutes exceeded 
the stress limits. For those cases considering growth only, the stresses 
near and on the inside surface of the chamber were highly compressive. 

However, when creep was included in the analysis, these same stresses were 
highly tensile, again illustrating the need for good creep data. 
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VI. PYROLYTIC GRAPHITE PRODUCT EVALUATION 


The quality of pyrolytic graphite which can be produced at the 
current stage of development was evaluated by a vendor survey and by labora¬ 
tory testing of free standing PG tubes made by a variety of techniques. Those 
aspects of PG quality of greatest interest for free standing PG thrust chambers 
are reproducibility and high internal pressurization strength. 

A. Vendor Survey 


The quality of PG depends on many factors such as furnace temperature 
and pressure, gas flow rate, mandrel material, and gas dilution by second¬ 
ary gases such as argon or hydrogen. A survey was made of the deposition 
techniques presently being used by the manufacturers of PG. Those vendors 
interviewed were as follows: 

1. Super Temp Corporation 
Santa Fe Springs, California 

2. General Electric Company 
Metallurgical Products Dept. 

Detroit, Michigan 

3. Pyrogenics, Incorporated 
Woodside, New York 

L. High Temperature Materials, Incorporated 
Lowell, Massachusetts 

5. Raytheon Company 
Research Division 
Waltham, Massachusetts 

Some of the techniques used by the manufacturers are considered by 
them to be proprietary and are not described in thi6 report. A summary of 
the results of the survey is given in Table X. Much of the development of 
PG has been for applications which have different requirements than those 
of free standing PG rocket thrust chambers. In particular, the requirement 
of high internal pressurization strength does not appear to have been a re¬ 
quirement for other applications, and hence the effect of the many deposition 
conditions on the rupture pressure of free standing PG had not been evaluated. 
The large number of possible combinations of the various deposition para¬ 
meters precluded a complete evaluation of all of the possibilities. However, 
for this program, the vendor’s judgment of those techniques which might most 
likely achieve high rupture pressure was solicited during the procurement of 
PG tubes for laboratory testing. Following the vendor survey, it was possible 
to evaluate the potential advantages of the deposition techniques suggested 
by the vendors, thereby selecting those techniques which appeared to have 
the greatest probability of success. 

B. Evaluation of Pyrolytic Graphite Tubes 

Two-inch diameter tubes of free standing pyrolytic graphite were 
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purchased from five vendors as part of a program to study the effects of 
various deposition techniques on rupture pressure. Since the primary goal 
was to maximize the rupture pressure, no limitations were placed on wall 
thickness. The vendors were asked to suggest any new or unusual deposition 
techniques which they thought might increase the rupture pressure. In addi¬ 
tion, FG tubes made by a variety of proven techniques were obtained. All 
tubes were 4 inches long, and they had a nominal inside diameter of 2 inches 
and a nominal wall thickness of 0.060 inch, although a number of tubes of 
wall thickness greater than 0.060 inch were produced. 

Since the object of the tube study was to investigate the rupture 
pressure of PG made by a variety of techniques, PG tubes produced by conven¬ 
tional techniques were not ordered from each of the vendors. Therefore, the 
strength of the tubes which were procured did not necessarily provide any 
comparison of the different vendors' overall capabilities. 

The tubes were pressurized to rupture with the test setup shown 
in Figure 91* The appearance of a typical tube after rupture is shown in 
Figure 92. The pressurizing medium was soluble oil and water. No axial 
load was applied to the tube during the rupture test. A number of tubes were 
instrumented with circumferential strain gages on the outside surface to 
measure the modulus of elasticity during pressurization and the residual 
stress after rupture and strain gage isolation. 

After rupture, 10X photographs of the surfaces and 50X photomicro¬ 
graphs of the c-plane edge were taken. The boron content of the boron-doped 
tubes were measured by the vendors. 

A description of the tubes which were purchased is given in 
Table XI. A discussion of the tubes provided by each vendor follows. 

1. Super Temp Corporation 

Four furnace runs were made by the Super Temp Corporation with 
eight mandrels loaded per run. The furnace conditions were kept uniform at 
a temperature of 4000°F and a pressure of 8.5 nim Hg. The variables studied 
were mandrel materials and mandrel machining. Only female mandrels were 
used, and the inside surface was machined with various lathe speeds and 
depths of final cut, as given in Table XII. All of the FG tubes produced 
by Super Temp were substrate nucleated. 

a. Effect of Mandrel Material 


It was found that the grain size of FG is related to the 
grain size of the graphite mandrel material. This i6 6hown by comparing 
Figures 93 through 96. Each of these figures show 10X photographs of the 
inside and outside surfaces and a 50X photomicrograph of a FG tube sample 
made on different mandrel materials. 

Tube S-103, 6hown in Figure 93, was deposited on a CS graphite 
mandrel (Maximum particle size = 0.030 inch). The grain size of Tube S-103 
was about 0.020 inch. However, superimposed on the grain matrix were a very 
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large number of nodules (anomalously large grains) up to 0.C40 inch in dia¬ 
meter on the inside surface. Tube S-103 ruptured at 370 psig. 

Tube S-102, shown in Figure 94, was deposited on an ATJ 
graphite mandrel (Maximum particle size = 0.006 inch). The grain size of 
Tube S-102 was about 0.015 inch. There were also a large number of nodules 
up to 0.060 inch. This tube ruptured at 380 psig, with a crack adjacent to 
an 0.060 inch nodule. 

Tube S-101, shown in Figure 95> was deposited on a Speer 
3^996 graphite mandrel (Maximum particle size = 0.003 inch). The grain size 
of the PG was about 0.015 inch. 

Tube S-104, shown in Figure $ 6 , was deposited on a Poco 
EP 192C graphite mandrel (Maximum particle size = 0.001 inch). The grain 
size of the PG was about 0.008 inch. 

The yield (percentage of uncracked tubes) and the rupture 
pressures of tubes deposited on CS and HIM mandrels were lower than for 
tubes deposited on ATJ, Speer, or Poco mandrels. One out o** two tubes de¬ 
posited were cracked when using CS and HIM mandrels, and the rupture pres¬ 
sures of the uncracked tubes were under 400 psig. 

Tubes deposited on ATJ mandrels were of intermediate 
strength (485 psig average rupture pressure) and the yield was four tubes 
out of four mandrels loaded. 

The average rupture pressure of tubes deposited cn Speer 
3499S graphite mandrels was 525 psig. For tubes deposited on Poco EP192C 
graphite mandrels, the average rupture pressure was about 590 psig, with 100 
percent yield in both cases. 

The results of these tests show that the attainment of 
high rupture pressure for thrust chambers which are too large for ATJ graphite 
mandrels may be difficult. A later furnace run was made by Super Tem? with 
infiltrated CS and HIM graphite mandrels. No improvement in production 
yield or tube rupture pressure was obtained. 

b. Effect of Nodule Size 


Earlier rupture tests of PG tubes, reported in Reference 
22 , had shown a correlation between small nodule size and high rupture pres¬ 
sure. A weak; correlation of the same type can be found in the rupture test 
results of the Super Temp tubes. One reason for this weak correlation may 
be that the nodules in these tubes usually occurred in clusters or strings, 
and the measurement of the nodule effect in terms of the maximum nodule diameter 
is probably not very accurate. 

It was hoped that the finer particle size and higher purity 
of Poco graphite mandrels might reduce the occurrence or the size of the 
nodules. This was not achieved during the four furnace runs by Super Temp. 
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However, visual inspection of the tubes deposited on Poco mandrels showed 
that the matrix of the PG tubes was fine grained and contained only small 
nodules. It still seems likely that improved PG strength should be obtained 
by using Poco mandrels. On the other hand, Poco graphite is not currently 
available in large enough sizes except for small thrust chambers. Several 
samples of PG with large nodules were polished on the c-plane edge to show 
the nodules and the point of their origination in cross section. By suc¬ 
cessive trials of polishing and examination under a microscope, it was 
concluded that probably all of the nodules in PG deposited on Poco graphite 
were originated at the graphite substrate. The cause of these nodules i6 
still unknown. However, it is believed that the cause is probably small 
particles of some sort which are on the surface of the mandrel and act as 
a nucleating source for the nodule. Strict mandrel cleaning procedures 
did not eliminate this problem, and it is probable that the dirt particles 
accumulate on the mandrel surface while the mandrels are being loaded into 
the furnace. The carbon black insulation used in the furnace is one possible 
source of the impurity. 

2. General Electric Company 

Two furnace runs were made by the General Electric Company 
with eight mandrels per run. The first furnace run was made at 3800°F to 
produce boron-doped PG with a boron content of about 0.5 percent. Six good 
tubes were delivered from this run. The second furnace run was made with 
deposition conditions (330C°F) which were expected to produce continuously 
nucleated PG. Five good tubes were delivered. All mandrels were ATJ 
graphite, polished with 600 grit paper. All deposition conditions except 
temperature are considered proprietary by General Electric Company. Data 
describing the General Electric tubes are given in Table XIII. 

The average rupture pressure of the boron doped tubes was 
85O psig, and the average for the continuously nucleated tubes was 1140 psig, 
including one tube (G203) which ruptured at 2100 psig, the highest rupture 
pressure so far achieved with a PG or PG alloy tube. 

a. Boron-Doped Pyrolytic Graphite 

Five tubes containing about 0.5 percent boron ruptured 
at an average rupture pressure of 850 psig. The grain size of these tubes 
was about 0.015 inch. The maximum nodule diameter in the uncracked tubes 
was 0.03 inch to 0.04 inch, while in the two cracked tubes it was 0.05 and 
0.10 inch. The high rupture pressure of these tubes was partly due to the 
thicker walls (about 0.080 inch) compared with the 0.060-inch wall thick¬ 
ness of most other tubes studied in this program. A correlation of maxi¬ 
mum nodule size with rupture pressure could not be established for the boron- 
doped tubes because four out of five tubes cracked near the same pressure, 

900 psig. 

The modulus of elasticity (about 6 x lO^psi) was higher 
than for most pyrolytic graphite. The residual stresses on the outside sur¬ 
face of six boron-doped tubes are given in Table XIII. The circumferential 
residual stresses on the outside surface ranged from -6530 to -8640 psi 
(compression). However, the measured circumferential residual stress on 


UNCLASSIFIED 

-48- 




AFRPL-TF-66-95 


UNCLASSIFIED 


Report 6l06 


the inside surface of one tube (G104) was only +4220 psi(tension), which is 
rather low, considering the wall thickness and the modulus of elasticity. 

The microstructure of Tube G102 shown in Figure 97 is 
typical of the boron alloy tubes. 

b. Continuously Nucleated Pyrolytic Graphite 

Five tubes of continuously nucleated PG were ruptured at 
pressures ranging from 690 to 2100 psig. The grain size of these tubes was 
about 0.010 inch, less than the grain size of the boron-doped PG tubes, 
although ATJ mandrels were used for both types of tubes. The grains also 
appeared to be harder, with more sharply defined grain boundaries, than the 
boron-doped PG. This can be seen by comparing the inside surfaces of Tubes 
G202 and G203 (shown in Figures 98 and 99) with the inside surface of boron 
PG Tube G102 shown in Figure 97* Photomicrographs of both ends of Tubes 
G202 and G203 are shown to illustrate the change in microstructure over a 
4-inch distance. One of the possible problems with continuously nucleated 
microstructure is the difficulty in maintaining the same microstructure over 
any length, since careful control of the local deposition conditions is 
required. It can be seen that there was a significant change in micro¬ 
structure from the top to bottom of Tubes G202 and G203, and in part, these 
tubes appeared to be substrate nucleated rather than continuously nucleated. 

The maximum nodule size in the continuously nucleated PG 
was 0.03 inch to 0.04 inch, and the absence of large nodules sometimes found 
in the boron PG tubes suggests the possibility that the growth of nodules 
is limited by the growth of new cones within the continuously nucleated PG. 

The modulus of elasticity in the a-direction was somewhat 
higher (about 5*0 x 10 & to 6.4 x 10 ^ psi) for the continuously nucleated 
PG than for substrate nucleated PG (about 4.5 x 10 ^). There was some scatter 
in the data, as shown in Table XIII. Whether or not the properties of con¬ 
tinuously nucleated PG are more susceptible to inconsistency than substrate 
nucleated PG is not known, although the scatter in the strengths of the 
4-inch long tubes seemed to indicate this. 

3 . Pyrogenics, Incorporated 

Three furnace runs were made by Pyrogenics, Incorporated. One 
tube was deposited in the first furnace run on a female mandrel, and contained 
controlled delaminations in a 0.l60 inch thick wall. This tube, designated 
P101, contained many large nodules, up to 0.10 inch in diameter, and the 
delaminations were quite nonsymmetrical. This tube was not ruptured because 
the ID was too large for the rupture test fixture. 

Tube P301 was a thicker wall (0.225 inch) tube deposited on a 
female mandrel in the second furnace run. There were many large nodules, up 
to 0.150 inch in diameter. The delaminations were nonsymmetrical and the ID 
was too small for the pressure test fixture. 
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Tubes P401 and P402 were deposited on a male mandrel. The 
wall of Tube P401 was about 0.30-inch thick and contained about six non- 
symmetrical delaminations. The maximum nodule diameter was 0.140 inch, on 
the outside surface. The inner layer of this tube ruptured at 100 psig, caus¬ 
ing leakage of the pressurizing liquid through the ends of the tubes. The 
microstructure and outside surface of Tube P401 are shown in Figure 100. 

The wall of Tube P402 was 0.120-inch thick and it contained two fairly uniform 
delaminations. The maximum nodule diameter on the outer surface was O .060 
inch. The inner layer of this tube ruptured at 40 psig. 

4. High Temperature Materials, Incorporated 
a. Controlled Delaminations 

Two furnace runs were made by High Temperature Materials, 
Incorporated, (HTM). In the first furnace run, three tubes with controlled 
delaminations were produced. The wall thickness was l/8 inch and the wall 
was divided into three approximately equal thickness cylinders by the two 
de laminations. The tubes were pressurized to rupture in the test setup 
(Figure 91 ) which applied hoop stress but no axial stress. Tube H101 rup¬ 
tured at 85 O psig, and Tube H102 ruptured at 1100 psig. Tube H103 was not 
ruptured because the ID was too small for the test fixture. 

The variation in strain on the outside surface with changes 
in internal pressure is shown in Figure 101 for Tube H101 and in Figure 102 
for Tube H102. It is shown that the outer shell does not start to carry much 
hoop tension until about 500 psig internal pressure has been applied. This 
is further illustrated in Figures 103 and 104, which show the variation of 
indicated hoop stress (Ec aa ) on the outer surface (from strain gage measure¬ 
ments) with theoretical hoop stress (Pr/t) for a solid, thin walled cylinder. 

At rupture. Tube H101 was carrying only 700 psi tensile 
stress on the outer shell, as contrasted to an average theoretical hoop stress 
for the entire tube of 6400 psi. 

The outer shell of Tube H102 carried a stress of 2200 psi 
at rupture, compared to an average theoretical hoop stress at rupture of 
85 OO psi. It is interesting to note that the curve of indicated hoop stress 
versus internal pressure is almost duplicated for both tubes, although one of 
the tubes (HlOl) ruptured at a lower pressure than the other. 

It appears that the controlled delamination concept is 
useful for extending the pressurization strength of PG chambers. Further im¬ 
provement is possible by making thicker walls (with more delaminationb) and 
perhaps by achieving a tighter fit between the lamina.. However, some pro¬ 
blems may arise when chamber shapes are deposited. The concept requires 
further evaluation and testing. 

bi Fibre-Reinforced Pyrolytic Graphite 

In the second furnace run, three PG tubes were deposited 
by HTM on a carbon fabric mandrel. The wall thicknesses were about l /8 
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inch. The rupture pressures were 600 psig for Tube H202 and 58 O psig for 
Tube H203. The third tube was not ruptured because the ID was too small for 
the test fixture. 

A photograph of Tube H202 is shown in Figure 105 and the 
microstructure and outside surface are shown in Figure 106 . The surface was 
quite uneven, compared to conventional PG tubes. Although the rupture 
strength was fairly good, the fibre-reinforced PG concept would require ad¬ 
ditional development. 


A number of furnace runs were made by the Raytheon Company to 
produce PG tubes of several different types. Two tubes were made in each 
run. A description of the Raytheon tubes is given in Table XIV. 

Two tubes were deposited in the first furnace run at a de¬ 
position rate of 17.2 mils/hr. Both tubes (R101 and R102) ruptured near 
400 psig. The wall thickness was about 0.060 inch. The microstructure of 
Tube R101 is shown in Figure 107 • 

Six tubes were deposited in a series of furnace runs at 
20 mils/hr. None of the six tubes had very high rupture pressure, ranging 
from 60 to 460 psig. 

Two tubes were deposited at 15 mils/hr to a wall thickness of 
about 0.090 inch. Tube R301 ruptured at 460 psig and Tube R302 ruptured 
at 5*+0 psig. 

Two tubes (R401 and R402), which were deposited at 30 mils/hr 
to a wall thickness of 0.090 inch, were ruptured at 230 and 380 psig, res¬ 
pectively. 

Four tubes containing about 1.5$ boron were produced. Two 
tubes were cracked, and two others ruptured at 320 and 380 psig. The micro¬ 
structure of Tube R502, which ruptured at 380 psig, is shown in Figure 107• 

Two tubes were deposited under conditions which were expected 
to produce high boron content ( 5 $ to 10$). However, the actual boron content 
achieved was only 1 to 2$. These tubes (R901 and R1001) were cracked and 
no rupture tests could be performed. The microstructure of Tube R1001 is 
shown in Figure 107. 

Finally, two tubes with a l/8-inch thick wall and containing 
one delamination were produced. The maximum nodule diameter in these tubes 
was 0.030 inch. 

Tube R1101 ruptured at 1260 psig. The stress on the outer 
surface at rupture was indicated by a strain gage to be 65 OO psi, compared 
to an average hoop stress, based on the entire thickness, of 11,100 psi. 

Tube R1102 ruptured at 1000 psig, at which time the stress 
carried by the outer surface was 58 OO psi, compared to an average hoop stress. 
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based on the entire thickness, of 8500 psi. 

The delaminated tubes produced by Raytheon were similar in 
strength to those produced by HIM, although they contained only one de- 
lamination instead of two. 

The rupture strength of most of the Raytheon tubes was low, 
with the exception of the two delaminated tubes. This is attributed to the 
unfavorable number and size of nodules. Raytheon is of the opinion that 
this may be due to a batch of ATJ graphite which was quite porous. In any 
case, no advantage could be established for any of the deposition parameters 
being evaluated. 

C• Conclusions 


1. The highest rupture pressures of any PG tubes tested were at¬ 
tained using a continuously nucleated microstructure. The application of this 
type of PG to thrust chamber fabrication will require additional process devel¬ 
opment to provide the carefully controlled microstructure throughout the 
chamber length and thickness. 

2. Pyrolytic graphite tubes containing carefully controlled de¬ 
laminations can currently be produced to carry about 1000 psig internal 
pressure at room temperature. This type of structure may encounter other 
problems, however, when tested in rocket firings. 

3. Boron alloys of PG as produced by the General Electric Company 
have shown higher tube rupture pressures than most substrate nucleated pyro¬ 
lytic graphites. 
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TABLE II 

THRUST FOR 1OO-P0UND THRUST CONFIGURATION 
«= 0*75 in. 

C* Efficiency - 9 % 

Shifting Equilibrium 


Propellant 

Combination 

Chamber 

Pressure 

P c 

Design 

Mixture 

Ratio 

o/f 

Sea Level 

Thrust 

e = 2 
e 

Space 

Thrust 

c * 40 
e 


(peia) 


(lbs) 

(lbs) 

N 2 0 u /0.5 N 2 H 4 - 0.5 UDMH 

50 

2.0 

19.5 

42.1 

100 

2.0 

51.9 

84.0 


300 

2.0 

181.0 

251.0 

F 2 /BA10l4 

50 

2.0 

20.9 

41.2 ! 


100 

2.0 

52.2 

82.1 


300 

2.0 

182.8 

244.0 ! 

F A 

50 

12.0 

21.0 

41.8 

100 

12.0 

52.6 

83.2 j 


300 

12.0 

l8l.O 

247.0 
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TABLE III 

WALL THICKNESSES OF 100-POUND THRUST CONFIGURATION 
Contraction Ratio = 7.14 


Thrust 

(lbs) 

Chamber 

Pressure 

(psia) 

Throat 

Diameter 

(in.) 

Chamber 

Radius 

(in.) 

100 

AH 

0.75 

1.00 


Chamber 


(in.) 


'3 

45 

0.06 
>8 0.08 


Throat 
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DIMENSIONS OF IOOO- AND 5000-POUND CONFIGURATIONS 

Shifting Equilibrium 

Expansion Area Ratio = 40 
95^ C* Efficiency 
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TABLE V 

WALL THICKNESSES OF 1000 -POUND THRUST CONFIGURATION 
Contraction Ratio = 4 


Thrust 

(lbs) 

Chamber 

Pressure 

(psia) 

Throat 
Diameter 
(inB.) 

Chamber 

Radius 

(ins.) 

Chamber 

Throat 


E9 

gjj | 

t 

(in.) 

1000 

50 

3.66 

3.6 

0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 

0.108 

0.162 

0.216 

0.288 

0.360 

0.432 

0.576 

0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 

0.055 

0.082 

0.101 

0.146 

0.183 

0.220 

0.293 

1000 

100 

2.6 

2.6 

0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 

0.078 

0.117 

0.156 

0.208 

0.260 

0.312 

0 . 4 i 6 

0.03 

0.045 

0.06 

0.08 

o.ro 

0.12 

0.16 

0.039 

0.058 

0.078 

0.104 

0.13 

0.156 

0.208 

1000 

300 

1.50 

1.50 

0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 

0.045 

0.067 

0.090 

0.120 

0.150 

0.180 

0.240 

0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 

0.022 

0.035 

0.045 

0.060 

0.075 

0.090 

0.120 
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Thrust 

(Its) 

Chamber 

Pressure 

(psia) 

Throat 

Diameter 

(ins.) 

Chamber 
Radius 
(ins.) 

5000 

50 

8.25 

8.25 

5000 

100 

5.85 

t 

1 

5.85 

5000 

300 

3.40 

3.40 


Chamber 


Throat 


0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.l6 


0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 


0.03 

0.045 

0.06 

0.08 

0.10 

0.12 

0.16 


0.247 


0.124 

0.186 

0.248 

0.330 

0 . 1>12 

0.495 

0.66o 


0.136 

O.170 

0.204 

0.272 
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TABLE VIII 

PARAMETERS FOR THE 5000-POUND THRUST CHAMBER SOAK BiEK ANALYSIS 

P ■ 100 psia 
c 

L* = 40 

Throat Diameter = 5*850 in. 


Injector 

Emissivity 

Emissivity 

Notes 

Chamber Wall 
Temperature 
Near Injector 

CF) 

Chamber 
Temperature 
Near Throat 

Cf) 

Film 

Cooling 

Notes 

Contraction 

Ratio 

Plotted 

in 

Figure 

0.25 

A 

5600 

5600 

F 

4 

55 

0.25 

A 

4500 

4500 

D 

4 

55 

0.25 

A 

2500 

4500 

E 

4 

55 

0.80 

C 

4500 

4500 

D 

4 

52 

0.25 

A 

4500 

4500 

D 

4 

52 

0.05 

B 

4500 

4500 

D 

4 

52 

0.25 

A 

4500 

4500 

D 

6 

54 

0.25 

A 

4500 

4500 

D 

4 

54 

0.25 

A 

4500 

4500 

D 

2 

54 


Notes : 

A. Emissivity of monel used in previous analysis 

B. Emissivity of polished pure silver = 0.02 to 0.03 
(Marks Mechanical Engineers Handbook 

McGrav Hill, New York 1958) 

C. Emissivity of injector covered with deposit 

D. Film cooling to inner‘Wall temperature of 4500°F throughout chamber 

E. Film cooling to inner wall temperature of 2500 °F for 5 inches from 
injector, and 4500°F elsewhere 

F. No film cooling 
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TABLE X (Continued 



Types of prerun treatments used or evaluated by several of 
the vendors include: flash coating, coating with colloidal 
graphite, high temperature mandrel soaking to eliminate out 
gassing, chlorine gas purification and others not specified 
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on vail thickness for Within 10 percent easily or 1 to 5 percent if required, 
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TABLE XI 


DESCRIPTION OF PYROLYTIC GRAPHITE TUBES 


Vendor 

No. of 
Tubes 

Description 

Super Temp Corp. 

32 

Eight per furnace run. Different mandrel 
materials and mandrel finishing proce¬ 
dures. Uniform furnace conditions. 

General Electric Co. 

8 

1 furnace run of boron-doped PG 


8 

1 furnace run of continuously nucleated 

PG 

Pyrogenics, Inc. 

1 

Controlled delaminations, female mandrel, 
0.l60-inch wall 


1 

Controlled delaminations, female mandrel, 
0.225-inch wall 


2 

Controlled delaminations, male mandrel 

High Temperature 
Materials, Inc. 

3 

Controlled de laminations, 0.125-inch wall 


3 

Fabric reinforced, 0.125-inch wall 

Raytheon Co. 

B 

17 mils/hr deposition rate, 0.06-inch 
thick 



20 mils/hr deposition rate, 0.06-inch 
thick 



15 mils/hr deposition rate, 0.09-inch 
thick 


D 

30 mils/hr deposition rate, 0.09-inch 
thick 


n 

1.5$ boron 


B 

High boron content 



Controlled delamination 


78 


A 
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SHIFTING EQUILIBRIUM 
EXPANSION AREA RATIO - 40 



1.4 1.8 2.2 2.6 3.0 3.4 


410 


390 


370 



OXIDIZER-FUEL RATIO 


FIGURE 3. Propellant Performance, F 2 /BA10l4 
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MEAN COEFFICIENT OF THERMAL EXPANSION - IO~° in./in 
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TEMPERATURE- °F 


FIGURE 11 . Thermal Conductivity in the a-Direction for Pyrolytic Graphite 

Analysis Standard 
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TEMPERATURE - °F 


FIGURE 12 . Thermal Conductivity in the c -Direction for Pyrolytic Graphite 

Analysis Standard 
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FIGURE lU» Modulus of Elasticity in the a-Direct ion for Pyrolytic Graphite 
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THROAT WALL TEMPERATURE 
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TIME - seconds 


FIGURE 17. Throat Wall Temperature vs. Time, 100-lh Thrust 
Configuration, NgO^/O.5 NgH^-0.5 UEMH, P c = 50 psia 
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FIGURE 19. Throat Wall Temperature vs. Time, 100-lb Thrust 
Configuration, NgO^/O.5 NgH^-0.5 UIMH, P c = 300 psia 
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FIGURE 23. Combustion Chamber Wall Temperature Distribution, 100-ib 
Thrust Configuration, NgO^/O.5 NgH^-0.5 UDMH, P c * 100 psia 
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OXIDIZER-FUEL RATIO 


FIGURE 25. Combustion Chamber Inside Wall Temperature vs. Oxidizer-Fuel 
Ratio, 100-lb Thrust Configuration, NgO^/O.5 NgH^-0.5 UEMH 
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TIME - seconds 


FIGURE 26. Throat Wall Temperature vs. Time, 100-lh Thrust Configuration, 

F 2 /BA10l4, P c = 50 psia 
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FIGURE 27. Throat Wall Temperature vs. Time, 100-lb Thrust Configuration, 

F 2 /BA10lU, P Q «= 100 psia 
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FIGURE 28. Throat Wall Temperature vs. Time, 100-lb Thrust Configuration, 

Fg/BA1014, P c = 300 psia 
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FIGURE 29. Combustion Chamber Wall Temperature vs. Time 
Configuration, Fg/BAlOlU, P c « 50 peia 
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FIGURE 31. Combustion Chamber Wall Temperature vs. Time, 100-lb Thrust 
Configuration, Fg/BAlOlU, P c = 300 psia 
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FIGURE 32. Combustion Chamber Wall Temperature Distribution, 100-lb 
Thrust Configuration, Fg/BAlOlU, P q = 100 psia 
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FIGURE 33 
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OXIDIZER-FUEL RATIO 


FIGURE 3^* Combustion Chamber Inside Wall Temperature vs. Oxidizer-Fuel 
Ratio, 100-lb Thrust Configuration, Fg/BAlOlU 
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FIGURE 35 • Throat Wall Temperature Comparison, 100-lb Thrust 
Configuration, Fg/BA10l4 and Fg/Hg 
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FIGURE 37. Combustion Chamber and Throat Inside Wall Temperatures 
vs. Oxidizer-Fuel Ratio, lOO-lb Thrust Configuration, Fg/Hg 
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FIGURE 39. Variation of Radial Thermal Conductivity with Throat Wall 
Temperature and Time, 100-lb Thrust Configuration, Fg/BAlOl^ 
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FIGURE 4l. Combustion Chamber Wall Temperature vs. Time, 1000-lb Thrust 
Configuration, NgO^/0.5 NgH^-0.5 UDMH, P c = 1-00 psia 
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FIGURE 45. Throat Wall Temperature vs. Time, lOOO-lb Thrust Configuration, 
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FIGURE I 4 .I 4 .. Combustion Chamber-Wall Temperature vs. Time, 5000 
Configuration, Fg/BAlOlU, P c = 50 psia 
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FIGURE 46. Combustion Chamber Wall Temperature vs. Time, 5000-lb Thrust 
Configuration, F 2 /BA10l4, P c = 300 psia 
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FIGURE 47. Throat Wall Temperature vs. Time, 5000-lh Thrust Configuration, 

Fg/BA10l4, P c = 100 psia 
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FIGURE 49. Dimensions of Thermal Model Chambers Used for Preliminary 

Heat Soak Back Analysis 
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FIGURE 54. Injector Soak Back Temperature vs. Contraction Ratio, 
5000- lb Thrust Configuration 
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FIGURE 55. Injector Soak Back Temperature vs. Combustion Chamber 
Inside Wall Temperature at Shutdown, 5000-lb Thrust Configuration 
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FIGURE 56 . Residual Stresses on the Surfaces of As-deposited 
PG Cylinders and Cones 
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FIGURE 57. Model for Axial Stress Analysis at Throat 
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FIGURE 58. PG Thrust Chamber Configuration Used for Residual Stress Measurements 
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FIGURE 6l. Chamber Operating Stresses vs. t/r, 100-lb Thrust Configuration, 

N 2 0 u / 0.5 NgH^ - 0.5 UBMH 
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FIGURE 63 . Throat Operating Stresses vs. Chamber Pressure, 

100-lb Thrust Configuration, NgO^/O.5 NgH^ - 0.5 UDMH 
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FIGURE 64. Chamber Operating Stresses vs. C* Efficiency, 
100-lb Thrust Configuration, Fg/BAIO^ 
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FIGURE 66. Chamber Operating Stresses vs. t/r, 
100-lb Thrust Configuration, F^/BAlOlU 
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FIGURE 72. Chamber Total and Incremental Circumferential Stresses vs. t/r, 
100 -lb Thrust Configuration, Fg/BA10l4 
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FIGURE 73« Chamber Total and Incremental Axial Stresses vs. t/r, 
100-lb Thrust Configuration, Fg/BAlOl^ 
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FIGURE 7k. Throat Total and Incremental Circumferential Stresses at 600 seconds 
vs. t/r,i100-lb Thrust Configuration, F^/BAlOl^ 
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FIGURE 7°» Throat Circumferential Operating Stresses vs. Firing Time 
100-lb Thrust Configuration, F o /BA10lU, t/r = 0.12 
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FIGURE 78 . Throat Circumferential Postrun Stresses vs. Firing Time 
100-lb Thrust Configuration, F o /BA10l4, t/r = 0.0 8 
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FIGURE 85 . Chamber Operating Stresses at 600 seconds vs. t/r, 
5000-lb Thrust Configuration, Fg/BA10l4 
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FIGURE 88. Operating Stresses, Controlled Wall Temperature Chamber, 
vs. t/r, 100-lb Thrust Configuration 
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FIGURE 91* Setup for Rupture Test of Pyrolytic Graphite 
















FIGURE 92. Typical Pyrolytic Graphite Tube After Rupture Test Failure 
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A. OD Surface (lOX) 



B. ID Surface (lOX) 



C. Cross Section (50X) 



FIGURE 93* Microstructure of PG Tube S-105, 

Rupture Pressure = 570 psig, National CS Graphite Mandrel 
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A. OD Surface (lOX) 




C. Cross Section (50X) 


FIGURE 9^» Microstrue Lure of Tube S-102, 

Rupture Pressure = 380 psig, National ATJ Graphite Mandrel 
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FIGURE 95. Microstructure of PG Tube S-101, 
Rupture Pressure = 680 psig, Speer 34995 Graphite Mandrel 
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A. OD Surface (lOX) 



B. ID Surface (lOX) 



C. Cross Section (}0X) 


FIGURE 96 . Microstructure of PG Tube S-10^ 

• Rupture Pressure = 900 psig, Poco Graphite Mandrel 
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B. ID Surface (lOX) 



C. Cross Section (50X) 


FIGURE 97* Microstructure of PG Tube G102 
Rupture Pressure = 9^0 psig, Boron-Doped 
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C. Top Cross Section ( 5 OX) 


D. Bottom Cross Section (50X) 


FIGURE 98 Micrestructure of PG Tube G202, Rupt Pressure = 69 O psia, 
Regenerative Microstructure, ATJ Graphite Mandrel 
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A. OD Surface (lOX) 


B. ID Surface (lOX) 


C. Top Cross Section ( 5 OX) D. Bottom Cross Section (50X) 


FIGURE 99* Microstructure of PG Tube G?03, Rupture Pressure = ?100 psia, 
Regenerative Microstructure, ATJ Graphite Mandrel 
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B, Cross Section (llX) 


FIGURE 100. Microstructure of PG Tube PUOl, 
Rupture Pressure = 100 psig, ATJ Graphite Mandrel 
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STRAIN, < - microinch / inch 


FIGURE 102. Strain on Outride Surface (from Strain Gages) vs. Internal Pressure, 

Tube H102 (Controlled Delaminations) 
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THEORETICAL HOOP STRESS - psl 


FIGURE 103. Indicated Hoop Stress on Outside Surface 
vs. Theoretical Hoop Stress, Tube H101 (Controlled Delaminations) 
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THEORETICAL HOOP STRESS - psi 


FIGURE IOU. Indicated Hoop Stress on Outside Surface 
vs. Theoretical Hoop Stress, Tube H102 (Controlled Delaminations) 
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FIGURE 105. HIM Fiber Reinforced Pyrolytic Graphite 
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A. OD Surface (lOX) 



B. Cross Section (llX) 


FIGURE 106. Microstructure of PG Tube H202, 
Rupture Pressure = 600 psig 
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